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Abstract 

In this thesis we investigate several aspects related to the theory of fiuctuations in 
the Cosmic Microwave Background. We develop a new algorithm to calculate the 
angular power spectrum of the anisotropics which is two orders of magnitude faster 
than the standard Boltzmann hierarchy approach (Chapter |^). The new algorithm 
will become essential when comparing the observational results of the next generation 
of CMB experiments with theoretical predictions. The parameter space of the models 
is so large that an exhaustive exploration to find the best fit model will only be feasible 
with this new type of algorithm. We also investigate the polarization properties of the 
CMB field. We develop a new formalism to describe the statistics of the polarization 
variables that takes into account their spin two nature (Chapter |^). In Chapter ^ 
we explore several physical effects that create distinct features in the polarization 
power spectrum. We study the signature of the reionization of the universe and a 
stochastic background of gravitational waves. We also describe how the polarization 
correlation functions can be used to test the causal structure of the universe. Finally in 
Chapter ^ we quantify the amount of information the next generation of satellites can 
obtain by measuring both temperature and polarization anisotropics. We calculate 
the expected error bars on the cosmological parameters for the specifications of the 
MAP and Planck satellite missions. 

Thesis Supervisor: Edmund Bertschinger 
Title: Professor 

Co-Supervisor: Uros Seljak 

Title: Harvard Smithsonian Center for Astrophysics Fellow 



2 



Acknowledgments 



The time has come to thank all the people that made this work possible. I would 
like to thank Ed for his support and for giving me the freedom to pursue my own 
interests and collaborate with other people. When I arrived at MIT three years ago 
he suggested I should talk to Uros. This proved to be an excellent suggestion. We 
started working on CMBFAST, not to enthusiastically at times, "well it is faster, so 
what?" we thought. Since then we have done several other things together and I 
think we had and are having a lot of fun. This is in my mind the most important 
thing. I have benefited immensely from talking and working with him, and hopefully 
this will continue into the future (if he has time to spend with the rest of us mortals 
now that he will become a distinguished professor). 

I want to thank the members of my thesis committee for their suggestions, in 
particular Ed and Uros for going carefully through the thesis to help improve my 
poor enghsh. At MIT I enjoyed the company other graduate students and postdocs, 
Greg, Jim, Lam and Rennan among them. 

I also want to thank David Spergel, several parts of this thesis were done in 
collaboration with him. It was a real pleasure to be able to work with him both 
during my visits to Princeton and over e-mail. I am very grateful for his support. I 
also want to thank John Bahcall for inviting me to Princeton for extended periods of 
time giving me the chance to work more closely with Dave and talk with the people 
at the IAS. 

Several people have made these years much more fun. My "roommates" Alejan- 
dro, Fernando, Lucas and Susana share among them extremely important qualities. 
Formidable cooking abilities, exquisite taste for movies, obsession for cleanness, pas- 
sion for sports and all other necessary characteristic of the perfect roommate. Patricio 
and our "ex-social chair" Pablo contributed to my well-being providing up to date 
party information and many fun weekday lunches. I also enjoyed keeping e-mail and 
phone contact with my old friends from the university and high school. Special thanks 
go to Fernando for his periodic updates on the absurd twists of argentinian politics 



3 



and his regular one line messages to learn how things were going, "JJ the great" for 
explaining me much more that I ever wanted to know about X-ray lasers and other 
perhaps more useful thoughts about life and Arturo, my friend from the "real world" 
with whom I have shared lots of things since our highschool days. 

My parents and brothers have been a constant source of support. I enjoyed greatly 
our regular phone calls and visits which remind me that there are people who care 
about me back home. Hopefully the large "spatial distance" that separates us will 
not get in our way. 



4 



Chapter 1 
Introduction 



Since the initial discovery of anisotropies in the temperature of the microwave back- 
ground radiation (CMB) by COBE in 1992 [|T|, the field has been transformed from 
a theoretical exercise to an active experimental area of research. The large amount 
of data from CMB experiments and other observations is beginning to revolutionize 
cosmology and has the potential of changing dramatically our view of the universe. 
In the next decade we might obtain answers to very fundamental questions such as 
the age, size and fate of the universe we live in. A few percent accuracy on most cos- 
mological parameters seems achievable in the very near future, when the next CMB 
satellite missions fly. Should the inflationary paradigm prove to be correct, detailed 
measurements of the anisotropies coming from a number of ground based, balloon 
borne and satellite experiments could be used to determine a large set of cosmologi- 
cal parameters such as the density of the universe (Qq) or the Hubble constant (Hq) 
with unprecedented accuracy 0, ^, ^. 

At the time the anisotropies in the microwave background were produced the 
universe was almost perfectly homogeneous. The fluctuations were at the 10~^ level 
and a linear approximation to the evolution equations suffices. This allows us to make 
very accurate predictions which are to be compared with experiments [§|, §], 0, §| • Both 
the physics and the equations to be solved are relatively simple. CMB anisotropies 
are very clean from the theoretical point of view. This distinguishes them from other 
more local probes of cosmology where nonlinear effects complicate the theoretical 
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Figure 1-1: Compilation of measurements of the CMB anisotropy power spectra 
(from Tegmark http://www.sns.ias.edu/~max/cmb/experiments.html). The theoret- 
ical prediction for COBE normalized SCDM is shown for comparison. 



interpretation of the data. 

A large number of other detections have followed the initial discovery by COBE 
g |Tg, 0, 0, ig 0, [TB|, 0, ig, [g, ^ |2T[. Figure |n] shows all the resuhs in 
the literature at the time this thesis was written (April 1998) in terms of the angular 
power spectrum of the temperature fluctuations, C/. The easiest way to understand 
the meaning of the angular power spectrum is to consider a small region of the sky 
(for example the simulation in Figure [2-11) . In that case the temperature map can 
be expanded in Fourier modes and the Ci is the power (mean square amplitude) as 
a function of the spatial frequency denoted with I. The rapid evolution of the field 
will make the figure outdated soon. We can expect new results from several ground 
and balloon borne experiments in the near future. The MAP satellite from NASA 
will be launched by the year 2000 and a few years later Planck from ESA should 
be launched too. The sensitivity and angular resolution of these satellites will be 
outstanding; MAP is expected to have a noise level of 20fiK per resolution element (of 
approximately 0.3° x 0.3°). This translates into an uncertainty for each multipole due 



to noise in the detectors of AQ ^ y 2/(2/ + 1) x {O.llnKf. The Planck satellite will 
be equipped with very sensitive bolometers in the highest frequency channels (above 
100 GHz). These channels will also have the best angular resolution, ~ 0.1°. For 
Planck the noise in the multipoles is expected to be AC/ ~ y^2/ (2/ + 1) x (O.Ol/xK)^. 
The high sensitivity of the experiments coupled to the accurate theoretical predictions 
are the primary motivation for using CMB anisotropies as a probe for cosmology. 

The impressive angular resolution and noise levels of the future satellite missions 
create certain technical problems for the data analysis. The first step in the analysis 
is to go from the time ordered data (temperature or temperature differences measured 
by the different detectors on the satellite at each time) to maps of the CMB sky. A 
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computationally feasible algorithm for doing this has recently been developed |22|. In 
the next step one wants to optimally measure the power spectrum of the anisotropics 
from the maps which have 10^ — 10^ pixels. Direct evaluation of the likelihood func- 
tion involves the inversion of a huge correlation matrix and is computationally too 
expensive. A practical algorithm solving this problem has recently been developed 
for MAP Hj. 

Once the power spectrum of the anisotropics has been measured we would like to 
find the model that best fits the data and the confidence intervals on the model pa- 
rameters. The power spectrum of the CMB anisotropics depends on a large number of 
parameters. We would like to learn about the total density of the universe Qq, Hubble 
constant Hq, baryon density Qb, the cosmological constant A and the amplitude and 
scale dependence of the initial fluctuations. The exploration of this vast parameter 
space that characterizes the models will require many calculations of the theoretical 
power spectrum. This is computationally too expensive with the usual Boltzmann 
codes which need to solve a large system of differential equations at a large number of 
wavenumbers. This computational barrier was overcome by the development of a new 
algorithm based on a line of sight integral p[ which we present in Chapter ^. This 
algorithm has proved very useful for the rapid exploration of the phenomenology of 
CMB anisotropics. The power spectrum for any cosmological model can be obtained 
in only a few minutes on a normal workstation. It has also been extended to include 
active sources (as opposed to the primeval potential fluctuations of standard models), 
so that it can be used to make accurate predictions of the anisotropics produced by 
topological defect models [P^|. The development of this new line of sight algorithm 
is one of the most important results presented in this thesis. 

There is still a need to develop ways of extracting additional information from 
the CMB maps not encoded in the temperature power spectrum, which could prove 
important if topological defects or other non-gaussian sources are shown to be the 
cause of the anisotropics. If inflationary models are correct, the power spectrum fully 
describes the statistical properties of the primary anisotropics, but secondary effects 
like gravitational lensing can leave interesting information in higher moments. These 
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effects tend to be small, at least on large angular scales. 

The CMB radiation field is also expected to be linearly polarized due to Thomson 
scattering of the photons with free electrons p5| , |2^, ^ 0. The polarization 



of the CMB can provide in principle the same information about our universe that will 
become available from the temperature maps. The main disadvantage of polarization 
is that the CMB is predicted to be only 1 — 10% polarized depending on the angular 
scale. Future satellite missions will have the required sensitivity and the most likely 
source of problems will be contamination from galactic foregrounds. The nature and 
level of the foregrounds for polarization remain uncertain and we may have to wait 
until the satellites fly for a definitive answer. There will also be several ground based 
experiments before the satellites, eg. They may help determine the level 

of foreground contamination to be expected for the satellite missions and may also 
detect cosmological polarization. 

The detection of polarization is challenging but nevertheless the rewards can be 
great. We have shown that polarization may be able to constrain the ionization 
history of the universe very accurately [^. We know that the light from the first 
generation of objects reionized the universe after recombination, but the time when 
this occurs remains very uncertain. The scatterings of the CMB photons with the 
electrons after reionization will leave a signature in the CMB polarization that could 
be our first direct detection of the epoch of reionization. We will also show that large 
angular scale correlations in the polarization can be used to test the causal structure 



of the universe and could provide a direct test of the inflationary paradigm [32 



Finally the pattern of the polarization vectors can be used to detect the presence of 



a stochastic background of gravitational waves |3^, 0. We will discuss all of these 
issues in Chapter ^j. 

The polarization of the CMB has the potential of providing very interesting tests 
for our theories of the universe. The study of the properties of the polarization is 
a unifying theme of this thesis. In Chapter |^ we introduce the physics of CMB 
anisotropy and polarization and we develop a statistical description for the CMB 
polarization field. Chapter |^ describes the algorithm for computing both temperature 
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and polarization spectra. In Chapter ^ we look at several ways in which polarization 
can provide extra information not available in the temperature power spectrum. In 
Chapter ^ we quantify the benefits of using polarization for the determination of 
cosmological parameters by the MAP and Planck satellites. 

The anisotropies in the CMB radiation field encode an enormous amount of infor- 
mation about the universe we live in. In the next years the main driver of the field 
will be the rapid progress of experimental research. The results of these experiments 
will be used in conjunction with the results of other observations to test our theories 
of the universe with an accuracy which was unimaginable just a decade ago. It is a 
great time for cosmology, the availability of large quantities of data is making all the 
difference and hopefully will provide several interesting surprises. 
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Chapter 2 

Statistical Treatment of CMB 
Fluctuations^ 



In this chapter we will develop the statistical description of the CMB radiation field. 
We will enphasize the spin 2 nature of polarization and the mathematical tools needed 
to analyze it. We will present a physical interpretation of the new variables introduced 
by making a connection with the theory of gravitational lensing and a brief summary 
of the mechanisms that can produce anisotropies and polarization in our universe. 

Figure |2-1| shows a simulated map of CMB anisotropies (both temperature and 
polarization) for the standard cold dark matter model (SCDM)^. The anisotropies 
can tell us about the state of the universe at a very early time, only 300.000 years 
after the big bang when hydrogen recombined and the universe became transparent 
to the CMB photons. Since then the photons have been able to travel an enormous 
distance, approximately 10^° light years. Thus the CMB anisotropies map the most 
distant parts of our universe at a very early stage. 

The main cause for the temperature anisotropies are fiuctuations in the photon 
density at recombination. We see more photons coming from what were denser re- 
gions, producing a hot spot in the map. Photons coming from different directions 

^Partially based on M. Zaldarriaga & U. Seljak, Phys. Rev. D 55, 1830 (1997). 

2 This model has flo = 1, = 0, Qb = 0.05, Ho = 50 km/sec/Mpc, a COBE normahzed 
power spectrum of initial fluctuations with a spectral index n — 1 and no stochastic background of 
gravitational waves. 
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also come from regions which were moving at different velocities relative to us; the 
consequent difference in the Doppler shifts they suffered also cause anisotropics. Grav- 
itational redshift is the last source of anisotropics. Fluctuations in the gravitational 
potential will lead to different amounts of gravitational redshift for photons coming 
from different parts of the universe. We can also consider the presence of a stochastic 
background of gravitational waves. They will also create temperature anisotropics as 
photons get redshifted or blueshifted by different amounts according to their direction 
of propagation relative to that of the gravity wave. 

Polarization is produced by Thomson scattering of the CMB photons and elec- 
trons. Thus it can only be produced when hydrogen is ionized. Furthermore only 
radiation that is anisotropic can lead to non zero polarization because by symmetry 
if the radiation incident on an electron is isotropic there is no preferred direction for 
the polarization of the scattered light. The sources of anisotropics described in the 
previous paragraph will create some polarization when Thomson scattering is present. 



2.1 Characterization of the Radiation Field 

The aim of this section is to develop the mathematical tools needed to describe the 
CMB anisotropics. We will emphasize the description of the polarization which is new 
to this thesis. The CMB anisotropy field is characterized by a 2 x 2 intensity tensor lij. 
For convinience we normalize this tensor so it represents the fluctuations in units of the 
mean intensity [lij = 61 /Iq), it is dimensionless. The intensity tensor is a function of 
direction on the sky n and two directions perpendicular to n that are used to define its 
components (§1,62). The Stokes parameters Q and U are defined as Q = (In — /22)/4 
and U = /12/2, while the temperature anisotropy is given by T = (In + /22)/4 
(the factor of 4 relates fluctuations in the intensity with those in the temperature, 
/ oc T^). In Figure we represent the polarization using vectors with magnitude 
P = ^/Q^~+T7^ that form an angle a = | arctan(t//(5) with ei. In principle the 
fourth Stokes parameter V that describes circular polarization would also be needed, 
but in cosmology it can be ignored because it cannot be generated through Thomson 
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scattering^. While the temperature is invariant under a right handed rotation in the 
plane perpendicular to direction n, Q and U transform under rotation by an angle ip 

as 



Q' = Qcos2iIj + Usm2ip 

U' = -Qsm2ij + Ucos2ij (2.1) 

where e'^^ = costp §1 + simp §2 and §2 = — simp §1 + costp 62- We can construct two 
quantities from the Stokes Q and U parameters that have a definite value of spin (see 
Appendix ^ for a review of spin- weighted functions and their properties), 

{Q±iUy{h) = e^^'^iQ±iU){h). (2.2) 



We may therefore expand each of the quantities in the appropriate spin-weighted 
basis 

Tin) = '^aT,imYim{n) 

Im 

{Q + iU){h) = ^ a2,im 2Yim{n) 

Im 

{Q-iU){h) = ^a-2,im -^Yimin)- (2.3) 

Im 

Q and U are defined at each direction h with respect to the spherical coordinate 
system {eg,e^). Using the first equation in ( |A.5D one can show that the expansion 



coefficients for the polarization variables satisfy 0^2,/™ = Ci2,i-m- For the temperature 
the relation is a^^^ = ax^i-m- 



■^The presence of a large magnetic field at recombination could generate a small component of 
circular polarization. 
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Figure 2-1: Simulated temperature and polarization map for SCDM (2.5" x 2.5" field). 
The polarization vectors are shown together with the map of the temperature. The 
temperature ranges from —IdQO/iK < T < ISlO/iK while the maximum amplitude 
of the polarization vectors is P = 12SijlK. 
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The main difficulty when computing the power spectrum of polarization in the 
past originated in the fact that the Stokes parameters are not invariant under rota- 
tions in the plane perpendicular to n (equation The usual approach starts by 
expanding all perturbations in Fourier modes with wavevector k. While Q and U 
are easily calculated in a coordinate system where the wavevector k is parallel to z, 
the superposition of the different modes is complicated by the behavior of Q and U 
under rotations. For each wavevector k and direction on the sky n one has to rotate 
the Q and U parameters from this k and n dependent basis where the calculation 
is done into a fixed basis on the sky before adding them. Only in the small scale 
limit is this process well defined. This approximation had always been assumed in 
previous work 0, 0, ^, ^ ^. However, one can use the spin raising and lowering 
operators 6 and d defined in Appendix to obtain spin zero quantities. These have 
the advantage of being rotationally invariant like the temperature and no ambiguities 
connected with the rotation of coordinate system arise. Acting twice with 6 , 6 on 
Q ±iU in equation (^73|) leads to 






nl/2 




E 



Im - 




(2.4) 



The expressions for the expansion coefficients are 




(/ + 2)! 



.-1/2 



/ 




dnY;jh)3\Q + tU)ih) 




(/ + 2)! 
(/-2)! 



1-1/2 




(2.5) 
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Instead of a±2,«m it is convenient to introduce their linear combinations 0, 



(lEM — — + '^-2,/m)/2 

a_B,«m = 'i{0'2,lm — a-2,/m)/2. (2.6) 



We can define two quantities in real space, 



l,m 

B{h) = Y.^B,imYim{n)- (2.7) 



Lm 



The temperature is a scalar quantity under a rotation of the coordinate system, 
T'{fi = TZfi) = T(n), where TZ is the rotation matrix. We denote with a prime the 
quantities in the transformed coordinate system. Equation (|2.5| ) shows that up to 
an /-dependent factor a±2,/m are the expansion coefficients of two spin zero quantities 
d 2(g + iU) and d ^{Q - iU). As a result E{h) and B (n) are also invariant under 
rotations. 

It is interesting to analyze the behavior of E and B under a parity transformation. 
We will consider the case where we reverse the sign of the x coordinate but leave the 
others unchanged. In spherical coordinates this amounts to changing the sign of cj) 
while 6 and r remain the same. Under this transformation e'^ = — and = so 
the Stokes parameters transform as Q'{n') = Q{n) and U'{h') = —U{fi). With the 



aid of equation (A. 2) we can show that 



d\Q + iUyin') = d\Q-iU){h) 

d^{Q-iU)\n) = d\Q + iU){fi). (2.8) 



This implies that (equations ^]5| p.6| and |2] 



E'{ft!) = E{h) 

B'{n) = -B{h). (2.9) 
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These two new variables behave differently under parity: while E remains unchanged 
B changes sign , in analogy with electric and magnetic fields. The sign convention 
in equation ( |2.6|) makes these expressions consistent with those defined previously in 
the small scale limit 0. 

To characterize the statistics of the CMB perturbations only four power spectra 
are needed, those for T, E, B and the cross correlation between T and E. The cross 
correlation between B and E or B and T vanishes because B has the opposite parity 
to T or E. The power spectra are defined as the rotationally invariant quantities 

CbI = 21 -\- I ^^'^*BM^B,lm) 



m 

1 



CI 



2/ + ^ ^ 



^Y.^O*TM<^E,lm) (2-10) 



in terms of which, 



/ * 


'm 




= CTlSl'l6m'm 




rn' 


'(^E,lm) 


= CElSl'lSm'm 




m' 




= CBlSl>l6m'm 


{^T,l' 


m' 




= CciSl'l6m'm 


{'^B,l' 


ml 


■Q'E,lm) 


~ \(^B,l'm'^T,lm 



The brackets (• • ■) are ensemble averages. 



f2.111 



2.1.1 Correlators in Real Space 

Sometimes it is useful to investigate the correlations in real space. These expressions 
are needed, for example, to write the likelihood function of a measured CMB map. 
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We rewrite equation (p^ ) as 



Im 

Qin) = -E 

t^(n) = — aB,lm ^i,;m(A) — iaE,lm ^2./m(fi)- 



(2.12) 



We have introduced Xi,/m(n) = ( 2>^/m+ -2'5^im)/2 and X2,im(n) = ( 2>^m- -2>^m)/2. 
They satisfy X^^^ = and X^^^ = -X2,/_m which together with aE,im = a*E^i-m 

and as^im = o*Bi-m ^lake Q and f/ real quantities. 



In fact, Xi imiiv) and X2,im(n) can be written in the form Xi^/m(n) = y (2/ + l)/47r 
i^i,/m(^) e^"^-^ and X2,/^(n) = ^(2/ + l)/47r F2,,^(^) e*'"^ where F(i,2),«m(^^) can be 



expressed in terms of Legendre polynomials 



(Z-2)!(/-m) 



m ! 



^(/ + 2)!(/ + 
sm 6 



sm 2 J 



{I - 2)\{l - m)\ m 



\ (/ + 2)!(/ + m)! sin^^ 
+ (/ + m)P;™i(cos6 



-(/ - l)cos^P;"(cos( 



(2.13) 



Note that F2,/m(6') = if m = 0, a necessary condition given that the Stokes param- 
eters are real. 

The correlation functions can be calculated using equations ( p.ll|) and (|2.12| ), 



(T(1)T(2)) = Y.CTiY.^VM)^yim{2) 

I '- m 

(g(i)g(2)) = 5:c^J^x*,„(i)Xi,,„(2) 

I '- rn 

{U{l)U{2)) = Y.Cei\£x;U^)X,,U2) 



Cbi 
Cbi 



E^2,/m(l)-^2,«m(2) 

m 

E^Mm(l)-^l,«m(2) 



■*A subroutine that calculates these functions is available 
http://arcturus.mit.edu/~matiasz/CMBFAST . 



at 
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(T(i)g(2)) = ECc. E oYL{i)XiMm 

I ^ m 

(T(l)f/(2)) = z^CcJE o>^L(l)X2,/™(2) 



(2.14) 



where 1 and 2 stand for the two directions in the sky ni and n2. These expressions 
can be further simphfied using the addition theorem for the spin harmonics, 



121 + 1 
An 



(2.15) 



where P is the angle between fii and n2, and (V'i,V'2) are the angles (eg, e^) at fii and 
n2 need to be rotated to become aligned with the great circle going through both 
points. In the case of the temperature, equation ( |2.15|) gives the usual relation. 



2/ -I- 1 



(2.16) 



For polarization the addition relations for Xi and X2^im, are calculated from 
equation ( |2.15| ), 



$:x*,^(i)Xi,,^(2) 

m 

E^2,Zm(l)-^2,«m(2) 
m 

EX*,^(1)X2,,„(2) 

m 
m 

E 0^m(l)-^2,«m(2) 



2/ + 1 

47r 
2/ + 1 
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21 + 1 
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21 + 1 



Att 
21 + 1 



Fi,i2{l3) cos 211)1 cos 2%l)2 - -F2,«2(/?) sin 2il)i sin 2il)2 
Fi,i2{l3) sin2?/;i sin2?/^2 - F2,i2{l3) cos2'?/'i cos2?/'2 
Fi,i2i/3) sin 2'ijji cos 2'?/'2 + -^2,/2(/?) cos 2'?/'i sin 2'ijj2 
Fi,o(/9) cos 27/^2 



— -i- 



Att 



Fi,;o(/3) sin 2^-2 . 



(2.17) 



We can equivalently write Fi,/o(/?) = ^Jil - 2)!/(/ + 2)! P^i/3). 

The correlations in equation ( p.l4|) with the sums given by equation (|2.17|) are 
all that is needed to analyze any given experiment. These relations are simple to 
understand the natural coordinate system to express the correlations is one in 
which the ei vector for each point (1 and 2) is aligned with the great circle connecting 
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the two directions, the e2 vectors perpendicular to the ei. With this choice of reference 
frame we have fSH, 



2/ -I- 1 

{Qr{mm) = E 



I 47r 



I 47r 



CbiF.MP) - CE1F2MP) 
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(T(i)g.(2)) = -j:^cciF,,om 

(T(l)[/,(2)) = 0. (2.18) 



The subscript r here indicates that the Stokes parameters are measured in this par- 
ticular coordinate system. We can use the transformation laws in equation to 
write {Q, U) in terms of {Qr, Ur). When we use equation (|2.18|) for their correlations 
we recover our final result given by equations (|2.14] ) and (p.l7|) . 



2.1.2 Small Scale Limit 

There are several temperature and polarization experiments being planned or built. 
Many of them will attempt to measure the anisotropics in small patches of the sky. 
In this limit the sky can be approximated as flat and spherical harmonics can be 
replaced by Fourier modes. This has important practical benefits because it allows 
the use of the fast Fourier transforms to reduce computational time. In this section 
we will study how polarization is described in this limit. In the small scale limit one 
considers only directions in the sky n which are close to z, in which case instead 
of spherical decomposition one may use a plane wave expansion. For temperature 
anisotropics we replace 

Y^aT^yUn) — / dH T(Z)e*^-^, (2.19) 



Im 



SO that 



T(n) = (27r)-2y dH T{l)e'^-^. (2.20) 
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The Fourier coefficients satisfy, 



{T{l)T{l')) = {27rf Cti S^'il-l'). 



(2.21) 



To expand s = ±2 weighted functions we use 



2^m 



2^m 




(27r)-2l^V^-^. 
r 



(2.22) 



They lead to the following expression 



{Q + iU){h) = -{2nfJdH [E{l)+iB{l)]-d^e 



{Q-iU){h) = -{2nf I dH [E{l)-iB{l)]^d^e'^-^ 



(2.23) 



From equation ([A.2|) we obtain in the small scale limit 



(2.24) 



where {l^ + ily) = /e*"^'. 

The above expression was derived in the spherical basis where ei = eg and 62 = e^, 
but in the small scale limit one can define a fixed basis in the sky perpendicular to 
z, e[ = e..j. and €.'2 = By. The Stokes parameters in the two coordinate systems are 
related by 



(g + tuy 
{Q - tuy 



e-'"^iQ + tU) 



(2.25) 
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Combining equations ( p.23| )- (|2^ ) we find 



Q'{e) = {2ny^ J dH [E{l)cosi2(Pi) - B{l)sm{2(l)i)]e'^-^ 

U'{0) = {2ny^ J dH [E{l)sm{2(f)i) + B{l)cos{2(f)i)]e'^-^. (2.26) 

These relations agree with those given in which were derived in the small scale 
approximation. As already shown there, power spectra and correlation functions for 
Q and U used in previous work on the subject [0, 0, ^ ^, Q can be simply derived 
from these expressions. 

In the small scale limit the correlation functions are (in their natural coordinate 
system denoted with an r in equation PTTH) , 



Ct{9) = J ^ e'^'^^^^' Cti 

CQi9) = /^e^'^-^^MC^,cos2(20^)+CBzsin2(20^)] 

Cu{e) = J ^ e^'^™^<^' [Cei sin2(20^) + Cbi cos2(20^)] 

Cc{0) = /^e^''^°^^'C7c,,cos(20^), (2.27) 



or equivalently, 

Ct{9) = p-^CTiMW) 

CQie) + Cu{e) = J ^-^ {Cei + Cbi) Mm 

Cq{9)-Cu{9) = J^^(Cei-Cbi)MW) 

Cc{9) = -J^^CciJ2m, (2.28) 

Ju are cylindrical Bessel functions. 

These relations can be inverted to obtain the power spectra in terms of the corre- 
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lations, 



Cti 


= 271 [ 




Jo 


Cei 


= 2n f 




Jo 


Cbi 


= 2n f 




Jo 


Cci 


= -2n 



271 / ede Ct{9) Jo{w) 



(2.29) 



2.1.3 Analysis of All-Sky Maps 

In this section we discuss issues related to simulating and analyzing all-sky polar- 
ization and temperature maps. This should be specially useful for future satellite 
missions 



10| , which will measure temperature anisotropies and polarization over 
the whole sky with high angular resolution. Such an all-sky analysis will be of partic- 
ular importance if reionization and tensor fluctuations are important, in which case 
polarization will have useful information on large angular scales (Chapter where 
Fourier analysis (i.e. division of the sky into locally fiat patches) is not possible. In 
addition, it is important to know how to simulate an all-sky map which preserves 
proper correlations between neighboring patches of the sky and with which small 
scale analysis can be tested for possible biases. 

To make an all-sky map we need to generate the multipole moments ot,;™,, o.E,im 
and aB^im- This can be done by a generalization of the method given in 0. For each / 
one diagonalizes the correlation matrix Mn = Ctu ^22 = Cei, ^12 = ^21 = Cci and 
generates from a normalized gaussian distribution two pairs of random numbers (for 
real and imaginary components of ai±m)- Each pair is multiplied with the square root 
of eigenvalues of M and rotated back to the original frame. This gives a realization of 
0'T,i±m and aE,i±m with correct cross-correlation properties. For aB,i±m the procedure 
is simpler, because it does not cross-correlate with either T or E, so a pair of gaussian 

1 /2 

random variables is multiplied with C^i to make a realization of aB,i±m- 

Once aE,im and aB^m are generated we can form their linear combinations a2.im 
and a_2.im- Finally, to make a map of Q{n) and f/(n) in the sky we perform the 
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sum in equation (|2.3|) . To reconstruct the polarization power spectrum from a map 
of Q{n) and U{fi) one first combines them in Q + iU and Q — iU to obtain spin 
±2 quantities. Performing the integral over ±2^/m (equation p.5|) projects out ±2saim, 
from which aE,im and as^^m can be obtained. It is important to remember that in this 
treatment Q and U are defined using the spherical coordinate system. 

Once we have the multipole moments we can construct various power spectrum 
estimators and analyze their variances. In the case of full sky coverage one may gen- 
eralize the approach in |Tl| to estimate the variance in the power spectrum estimator 



in the presence of noise. We will assume that we are given a map of temperature and 
polarization with Npi^ pixels and that the noise is uncorrelated from pixel to pixel 
and also between T, Q and U. The rms noise in the temperature is ax and that in Q 
and U is ap. If temperature and polarization are obtained from the same experiment 
by adding and subtracting the intensities between two orthogonal polarizations then 
the rms noise in temperature and polarization are related by cr|i = ap/2. 

Under these conditions and using the orthogonality of the sY/m, we obtain the 
statistical property of noise, 

A 2 

// noise \* noise \ _ ___Tx X 

\\(^T,lm) °'T,Vm'l " divdmra' 

pix 

noise \* noise \ _ STTCTp ^ 
^^pix 

// noise \* noise \ _ ^ c 

U«-2,/mJ «-2,Z'm'/ " ^TT ('ii'fJmm' 

^^pix 

((«-tL)*«w) = 0- (2-30) 

By assumption there are no correlations between the noise in temperature and polar- 
ization. With these and equations ( p.6| ) and ( p. Ill ) we find 
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{^*B,l'm'^E,lm) — {(^*B,l'm'^T,lm) " 0- 



(2.31) 



For simplicity we characterized the beam smearing by e'^'^*/^ where cxb is the gaussian 
size of the beam and we defined w^p = Ana^ p/Npi^ f^, [Til • 



The estimator for the temperature power spectrum is [11 
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Tl 



E 



\0'T,lm\ 

21 + 1 



— Wq 



e b 



(2.32) 



Similarly for polarization and cross correlation the optimal estimators are given by, 



Cei 
Cbi 
Cci 



E\(^E,l'm\ —1 
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2/ + 1 
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21 + 1 



■''*E,lm'^T,lm + O'E.lmO'T,, 
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2(2/ + 1) 
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(2.33) 



The covariance matrix between the different estimators, Cov(XX') = ((X — 
{X)){X' — {X'))) is easily calculated using equation (|2.31|) . The diagonal terms 
are given by 



Cov(C'|,) 

Cov(C|,) 
Cov(C'S,) 



21 + 1 

^(C.. + Ve' 



iCBi + Wp'e' ^"f 



21 + 1 
1 

21 + 1 



Cci + {Cti + w^'e' ''^){Cei + Wp'e' 



(2.34) 



The non-zero off diagonal terms are 



CoviCriCEi) 
CoYiCriCci) 
CoviCEiCci) 



21 + 1 
2 

21 + 1 
2 

21 + 1 



-a 



Cl 



CciiCTi + Wr'e' 



Cci{CEi + Wp^e^ 



(2.35) 
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Note that the theoretical analysis is significantly more complicated if all four power 
spectrum estimators are used to deduce the underlying cosmological model. For 
example, to test the sensitivity of the spectrum to the underlying parameters one uses 



the Fisher information matrix approach |12, 03]. If only temperature information is 



given then for each I a derivative of the temperature spectrum with respect to the 
parameter under investigation is computed and this information is then summed over 
all / weighted by Cov^^(C'|,;). In the more general case discussed here instead of a 
single derivative we have a vector of four derivatives and the weighting is given by 
the inverse of the covariance matrix, 

«« = i:i:^Cov-'(C,,C,,)^. (2.36) 

I X,Y ^'^1- ^^3 

where aij is the Fisher information or curvature matrix, Cov^^ is the inverse of the 
covariance matrix, Si are the cosmological parameters one would like to estimate and 
X, Y stands for T, E, B, C. For each / one has to invert the covariance matrix and sum 
over X and Y , which makes the numerical evaluation of this expression significantly 
more involved. 



2.2 Understanding E and B Polarization 

In the previous sections we have discussed in detail the mathematical formalism 
needed to characterize the anisotropies in the CMB. In this subsection we want to 
take another look at the new polarization variables we introduced, E and B. By 
doing so we will develop further intuition as to what E and B physically mean. We 
will also study the power spectra and correlation functions for SCDM to point out the 
similarities and differences between the temperature and polarization anisotropies. 

2.2.1 E and B in Real Space 

The main characteristic of the polarization field is that it is a spin-2 field. The new 
variables E and B allow us to describe it more conveniently in terms of two spin zero 
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quantities. So far we have done this in Fourier space (more generally in /-space when 
not working in the small scale limit). In this section we want to explore E and S-type 
polarizations directly in real space. We will work in the small scale limit to make the 
notation simpler; the generalization to all sky modes is straightforward. 
As in equation ( ^.2UD we can define 



Eie) = {2n)-^ J dH e'^-^ E{1) 

B{e) = {2n)-^ J dH e'^-^ B{1). (2.37) 



These two quantities describe completely the polarization field. Figure |2-2| shows the 
E polarization field and the polarization vectors for the same simulation shown in 
Figure |2-1| . There is only E-tjpe polarization associated with SCDM because this is 
the only pattern that is produced by density perturbations and we did not include a 
stochastic background of gravity waves. We will discuss this point further in the next 
section and in Chapter ^. 

In Figure we can see that hot spots of the E map correspond to points with 
tangential polarization patterns (negative Qr)- We find radial polarization patterns 
around the cold spots of E. The polarization pattern in our simulation did not have 
any i?-type polarization so to get a better intuition as to what B-tjpe polarization 



means we can take the same E field in Figure |2-2| and pretend that it was actually B- 



type polarization. The polarization vectors are shown in Figure they correspond 
to the ones in Figure |2-2| rotated by 45°. Hot and cold spots of the B field correspond 
to places where the polarization vectors circulate around in opposite directions. It is 
clear from this Figure that the polarization pattern is not invariant under reflections. 
The distinction between E and -B-type of polarization is not in the size or orientation 
of the polarization vector at a given point but in the relative orientations of the 
vectors, namely the pattern. 

We can stress the characteristics of E and -B-type of polarization by looking at 



Figure p-4| which summarize what we found in the simulated maps. Places on the 



sky where the polarization pattern is tangential correspond to hot spots of E while a 
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radial pattern will produce a cold spot. The patterns that give rise to B are different 
in nature; while the radial and tangential patterns are symmetric under reffections 
those that create B are not. One pattern transforms into the other under a parity 
transformation. This is precisely what we mean when we say 5 is a pseudoscalar, if 
the pattern of polarization vectors on the sky was such that we had a hot spot of B 
when we looked at the mirror image of that pattern we would have a cold spot of B. 



We can use Figure |2-4| to argue why density perturbations cannot create 5-type 
polarization while gravity waves can, although the detailed mathematical proof is 
left for Chapter ^ We can consider only one Fourier mode of density perturbations 
because an arbitrary density field can always be expanded in Fourier modes. In 
the linear regime the total B generated by all the modes will be the sum of the 
contributions by each mode. We will prove the total -B-type polarization is zero 
because the B polarization generated every mode is zero. 

Figure p2-5| illustrates our point. In the upper half we consider the perturbations 
induced by a single density mode. The problem has two important symmetries, 
symmetry under rotation around k and symmetry under a reflection about any plane 
containing fc. It is clear that these two symmetries do not allow for a non zero U in 
the coordinate system shown in the figure. The polarization vectors must be either 
along eg or e^. The two possible orientations of the polarization consistent with the 
symmetries are shown in the upper right hand side. This can be shown to imply that 
5-type polarization is zero. We first focus on the particular case when the direction 



of observation coincides k because it is easier to visualize. We can use Figure \2-4\ with 
to recognize the polarization patterns that are allowed around k. The two -B-type 
polarization patterns cannot be produced because they are not invariant under parity. 
We leave the generalization of the argument to directions of observation different than 
k for later, when we introduce a way of calculating E and B directly in real space. 



31 



Figure 2-2: Simulated polarization map for SCDM(2.5'' x 2.5" field). The polarization 
vectors are shown together with the E-tjpe polarization. The E field varies between 
—liOliK < E < IbS/iK, the largest polarization vector has an amplitude P = 12SiiK. 
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Gravitational waves do not possess either rotational or reflection symmetry. We 
illustrate this in Figure |2-5| using the deformation of a ring of test particles in the 



X — y plane when a gravitational wave traveling along the z passes by. The influence 
of a gravitational wave does not have the symmetry under reflections that made zero 
the i?-type polarization for density perturbations and will thus produce some B. 
We can combine equation ( p. 37 ) with the deflnition of E{1) and B{1) in terms of 



Q and U, 

E{1) = J (fe \Q{e) cos(20j) + U{e) sin(20O] e~'^ ^ 

B{1) = J (fe [U{e)cos{2(j)i) -Q{e)sm{2(j)i)] e-'^-^, (2.38) 

to obtain an expression relating E and B in real space directly to Q and U, 

E{d) = - J (fe' uj{e) \Q{e') cos(20) - U{e') sin(20)] 

= - j (fe' uj{e) Qr{e') 

B{d) = - J Se' 00(6) [U{d') cos(20) + Q{d') sin(20)] 

= - / uj{e) Ur{e'). (2.39) 



The variables {6, 0) are the polar coordinates of the vector — 0' . In equation ( p.39|) 



Qr and Ur are the Stokes parameters in the polar coordinate system centered at 6. 
For example if is zero, Qr = cos 20' Q{0') — sin 20' U{6') and = cos 20' U{6') + 
sin20' Q{e'). The window can be shown to be u{e) = l/nO^ {9 ^ 0), uj{e) = {9 = 
0). We can read directly from equation ( p^.39 ) the relations between the pattern of 



polarization and the sign of E and B, for example a tangential pattern of polarization 
which has negative Qr will produce positive E (hot spot of E). 

From equation ( p.39[ ) we can understand how we achieve two rotationally invariant 
quantities out of the spin 2 fleld: to get E{0) and B{6) we average the values of Qr 
and Ur respectively, over circles centered at 0. Each circle is weighted by uj{6). By 
construction these two quantities are rotationally invariant: the Stokes parameters 
Qr and Ur do not depend on the coordinate system, they are deflned relative to the 
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Figure 2-3: Simulated polarization map for SCDM(2.5° x 2.5° field) where the E'-type 
polarization has been changed into S-type. The polarization vectors of the this figure 
are rotated 45° respect to the ones in previous ones. 

— 6' vector and the weight function uj is also rotationally invariant. We are giving 
the same weight to all the points in each circle and we are using the Stokes parameters 
defined in their natural coordinate system. The variable B is clearly a pseudoscalar 
because it is the average of Ur and Ur changes sign under parity. 

We can use Equation ( p.39|) to prove that i?-type polarization is zero in any 



direction for a single density mode. We had argued that this was the case when the 
direction of observation and k where the same. In general the two directions form an 
angle and k will intersect the plane perpendicular to n at a point different from n. 
We illustrate this in Figure |2-5| where we considered the case in which the polarization 
vectors are radial around fc, but the argument is the same if they were tangential (the 
only allowed patterns at a given distance away from k). To compute B we integrate 
Ur along circles centered at n, but the polarization pattern has reflection symmetry 
across the line connecting k and n which implies that on opposite sides of the circle 
Ur has different signs and the integral cancels. 

There is nothing special about the weight oj{9) as far as constructing a scalar and 
a pseudo-scalar, but this choice of weight has several important properties. It makes 
the correlation functions of E and B the closest to that of Q and U . For example 
with this choice Cq{6) + Cu{P) = Ce{6) + Cb{0)- Another property is that this choice 
of weight preserves the nature of white noise. White noise in Q and U becomes white 
noise in E and B. On the other hand this choice of window requires integration over 
all 9 and is non-local. 
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Another expression for E{6) and B{0) can be obtained in terms of derivatives, 

E{e) = V'\dl-dl)Q{e) + v-\2d^dy)u{e) 

B{e) = w-\dl-dl)u{e)~v-\2d,dy)Q{e). (2.40) 

This equation makes the point that the extra / factors that would be introduced in 
the power spectra by the derivatives are compensated by the inverse Laplacian. This 
is what we mean when we say that the spectrum of E and B are the closest to that 
of Q and U . The inverse Laplacian is responsible for the non-local nature of E{6) 
and B{6) evident also in equation (p.39|) . 



2.2.2 Analogy with Weak Lensing 

There is a very close analogy between the E and B formalism we are using to describe 
CMB polarization and the mathematical framework used to study weak lensing. As 
photons travel across the universe their trajectories are deflected. As a result the 
observed intensity is related to the true one by Iobs{G) = ItmeiG + ^G)- This mapping 
will induce ellipticities in the images of the observed galaxies. The deformation tensor 
^ij {hi = x,y m. the small scale limit) can be calculated in terms of a projected 
gravitational potential 0proiB 

= didj(j)proj. (2.41) 

The deformation tensor is usually decomposed into trace k and the two compo- 
nents of the shear 71 and 72, = — ^yy = —i^ + li, ^xy = ^yx = ~72- The 
trace n is also the dimensionless projected mass density. The shear part of the tensor 
can be obtained from the observed ellipticities of the galaxies . 



The two quantities 7(1,2) are the analogues of {Q, U) for the CMB polarization 
field. They also form spin two variables 71 ± ^72. Equations ( |2.4CI| ) and ( |2.41| ) can 



be used to show that the E variable in this case is nothing but the dimensionless 



^The interested reader may find detailed expressions for the projected potential and other quan- 
tities in ITl. 
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Figure 2-4: Polarization patterns that lead to positive and negative values of the E 
and B polarization fields. The Stokes parameters are measured in the polar coordinate 
system centered at the cross. All four patterns are invariant under rotation but the 
two patterns that generate B are not invariant under refiections. 



Figure 2-5: Upper half: A single mode of density perturbations has symmetries under 
rotation around the k axis and reflection about any plane containing it. Only I and Q 
Stokes parameters can be present in this reference frame and thus no B polarization is 
created. The integral along a circle of Ur, used to calculate B directly in real space, is 
zero because the pattern of polarization is symmetric under reflections across the line 
containing k and n. Gravitational waves do not have these symmetries as illustrated 
by the deformation suffered by a ring of test particles as wave traveling along z passes 
by. 

projected mass density n, 

E^^^O) = V-\dl - dl) 7i(0) + V-'my) 72(0) = VVpro,(0) = <e) (2.42) 

Thus in the case of weak lensing E has an important physical: interpretation it is the 
projected mass density. This is clear in Figure |2-2| : hot spots of the E fleld are the 
analog of high mass concentrations, so we see in the polarization vectors the typical 
orientation of the shear that we see in the deformation of galaxies around clusters. 
The analogous calculation of B for weak lensing is identically zero, as it is the z 
component of the curl of the gradient of the projected potential. 



2.3 Correlators for SCDM 

In this section we present the correlators for COBE normalized SCDM to point out 
some of the properties of CMB anisotropies and to understand their physical origin. 
Hydrogen recombination plays a crucial role in the generation of CMB anisotropies. 
As the universe expands it cools and around 3 x 10^ years after the big bang hydro- 
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Figure 2-6: Intensity patterns incident on an electron and the resulting polarization 
of the scattered light in the direction perpendicular to the page. The dashed line 
indicates a smaller intensity while the full lines represent an excess in the number 
of photons. A uniform intensity field or a dipole pattern will produce no polariza- 
tion. The rods around each diagram represent the direction of the polarization of 
the scattered light which come from that particular direction before the scattering. 
The total scattered light can be obtained by "adding" the four rods whose length is 
proportional to the amount of photons incident on the electron from each direction. 
Only the quadrupolar pattern will induce some polarization, because the summed 
lengths of the top and bottom rods differs from that of the left and right rods. 

gen atoms are able to recombine to form neutral hydrogen and the universe becames 
transparent to the CMB photons. Recombination occurs very quickly. Before recom- 
bination photons and electrons scatter very efficiently forming a single fluid, while 
after recombination the photons are free to travel to the observer. The photons we 
detect coming from a particular direction come from a very small region of our uni- 
verse (of size ~ 50 h^^ Mpc) at a distance D ~ 6000 Mpc away from us. The 
collection of all this regions in the universe is called the last scattering surface. 

In order to understand the main features in Figure P-l| one can work in the thin 
scattering surface approximation, where one assumes recombination occurs instanta- 
neously. In this approximation the final temperature fractional anisotropy in direction 
n on the sky is 

r(n) = ^1,, - n ■ vrI., - i fj drk.h^h^ (2.43) 

5r and vr stand for the relative density perturbation and velocity of the photon- 
baryon plasma while hij describes the perturbation in the metric. Conformal time 
is denoted with r; Tq and correspond to the times today and at recombination, 
respectively. The first two terms are evaluated at the last scattering surface and the 
third term is an integral along the line of sight. The anisotropics have three distinct 
origins: an overdensity of the photon baryon fiuid somewhere on the last scattering 
translates into an increase in the measured temperature when we observe in that 
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particular direction (first term). Furthermore, if tlie fluid is moving there will be an 
additional Doppler shift contribution (second term). Finally, gravitational redshift 
will change the photon temperature as photons travel toward us (third term). Note 
that this effect can take place anywhere along the path. 

In order to generate polarization we need Thomson scattering between photons 
and electrons, which means that polarization cannot be generated after recombination 
(if there is no reionization) . But Thomson scattering is not enough, the radiation 



incident on the electrons must also be anisotropic. Figure |2-6| illustrates this point. 
If the radiation incident on the electron is isotropic, by symmetry there can be no net 
polarization after the scattering. We consider a dipole anisotropy next: in the figure 
we illustrate the case where the intensity of the incident radiation is higher from the 
top and lower from the bottom, with the average intensity incident from the sides. 
The scattered radiation from photons incident from either the bottom or top will be 
polarized in the horizontal direction, while that coming from the sides will vertically 
polarized. The total degree of polarization is obtained by looking at all the scattered 
light, which in this case is unpolarized because the excess of radiation coming from 
above (relative to the sides) is compensated by the smaller amount coming from the 
bottom. Since a dipole pattern of anisotropies will not create any polarization, we 
need a quadrupolar pattern of intensity: an excess in photons must come from both 
top and bottom. 

In the context of CMB, velocity gradients in the photon-baryon fluid will create 
the quadrupole that generates polarization. The photons that are scattered from a 
given electron come from places where the fluid has velocity v. Because of the tight 
coupling between photons and electrons the photon distribution function has a dipole 
term Ti = fi- v. Furthermore, gradients in the velocity field across the mean free path 
of the photons (Ap) create a quadrupole T2 = Xpn^n^diVRj in the photon distribution 
as seen in the rest frame of the electron. The velocity of the fluid will only create a 
dipole; one needs at least a velocity gradient to create a quadrupole. 

The scattered radiation field is given by (Q+iU) = —3/4aT J dQ' /A'K{m-h')T2{h') oc 
XpTcn^mWiVjlr,, where ax is the Thomson scattering cross section and we have written 
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Figure 2-7: Two dimensional universe analogy. A particular density mode is oscillat- 
ing in time before decoupling. At decoupling the photons are free to escape and have 
traveled a distance D to reach the observer (at the center of the circle). Although 
the density field has continued to evolve after recombination, the photons carry the 
information of the state of the plasma at the time the last scattered. We would see a 
pattern of hot and cold regions in this one dimensional sky because photons arriving 
from different directions come from places with varying photon energy density. This 
pattern is shown as a dotted line on the right side of the circle. 

the scattering matrix as P(m, n') = — 3/4o"r|m ■ fi'p with m = ei + ie2 . In the last 
step we integrated over all directions of the incident photons n'. As photons decouple 
from the baryons their mean free path grows very rapidly, so a more careful analysis 
is needed to obtain the final polarization0|TI 



(g + iU){h) ^ O.nAnm'm^diVjlr, (2.44) 

where At^, is the width of the last scattering surface and is giving a measure of the 
distance photons travel between their last two scatterings. The appearance of m^m^ in 
equation ( |2.44|) assures that {Q + iU) transforms correctly under rotations of (ei, 62). 

We can combine equations ( |2.40| ) and ( |2.44| ) to understand why B-tjpe polar- 
ization is zero for anisotropics produced by density perturbations. We take n = z 
and m = X + iy. We then have V^-B oc {dy — d1)U + 2dxdyQ which gives V^B oc 
V^(z- V X v) and is zero because the velocity field produced by density perturbations 
is irrotational, v = Vtp and V x Vip = 0. On the other hand V^-E oc V'^{dxVx + dyVy). 

Before recombination Thomson scattering keeps the photons and baryons tightly 
coupled. They form a fluid with pressure provided by the photons and inertia by the 
baryons. This fluid supports the analog of acoustic oscillations where both the density 



and velocity are oscillating functions of time |T^, |T8[. The density is proportional to 
cos^Cskr) while the velocity to sin^Cskr), Cg is the sound speed. After hydrogen 
recombines the photons are free to travel to the observer without further scatterings. 



^The velocity in this equation is in the conformal gauge. 
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Figure 2-8: Power spectra for COBE normalized SCDM. The E and the T — E 
spectra have been rescaled for convenience. 



We illustrate the important ingredients in the generation of the anisotropies with 
a two dimensional analogy in Figure |2-7| . We show only one density mode which 
is oscillating as a function of time before recombination. Suddenly the hydrogen 
recombines when this mode was in a particular phase in the oscillation. The photons 
travel freely a distance D to the observer who sees more photons coming from the 
what were denser regions at recombination. The pattern of hot and cold temperature 
produced by this mode is not a perfect sinusoid because we are intersecting a plane 
wave with a circle, but it has a typical period. The contribution peaks at I ~ kD. 



In Figure |2-8| we show the three power spectra needed to characterize the CMB 
anisotropies. The oscillations in the CMB spectra can be analytically understood 
within this picture. Waves of different wavenumber k are at different phases in their 
oscillations at this particular time, which translates into an oscillating amplitude for 
the modes as a function of k. Physical size and angle in the sky are related by 
the angular diameter distance to recombination, so the oscillating amplitude of the 
modes also produces acoustic peaks in / space. On small scales (/ > 1000), when 
the wavelength of the perturbation becomes comparable to the mean free path of the 
photons prior to recombination, the anisotropies are suppressed by photon diffusion. 
Photons can diffuse out of density peaks, thereby erasing the anisotropies. 

The curves in Figure |2-8| illustrate the differences between temperature and po- 
larization anisotropies. The large angular scale polarization is greatly suppressed. 
Correlations over large angles can only be created by the long wavelength pertur- 
bations, but these cannot produce a large polarization signal because of the tight 
coupling between photons and electron prior to recombination. Multiple scatterings 
make the plasma very homogeneous; only wavelengths that are small enough to pro- 
duce anisotropies over the mean free path of the photons will give rise to a significant 
quadrupole in the temperature distribution, and thus to polarization. 
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Figure 2-9: Correlation functions in real space for COBE normalized SCDM model. 
The spectra have been smoothed with a 9fwhm = 0.2° corresponding to the beam size 
of MAP. 

On subdegree angular scales temperature anisotropy, polarization and their cross 
correlation show acoustic oscillations (Figure |2-8| ) , but in the polarization and cross 
correlation spectra the peaks are much sharper. The polarization is produced by 
velocity gradients of the photon-baryon fluid at the last scattering surface (equation 
|2.44|) . The temperature receives contributions from density and velocity perturbations 
and these partially cancel each other making the features in the temperature spectrum 
less sharp. The dominant contribution to the temperature comes from the oscillations 
in the density, which are out of phase with the velocity. This explains the difference 
in location between temperature and polarization peaks. The extra gradient in the 
polarization signal, equation ( p.44| ), explains why its overall amplitude peaks at a 
smaller angular scale. The fact that the polarization field has relatively more small 



scale power is evident when we compare the the T and E fields in Figures and 



2-2 



Figure |2-9| show the correlation functions in real space. The spectrum has been 
smoothed with a Ofwhm = 0.2° gaussian, similar to the MAP beam. An interesting 
point is that both polarization auto-correlation functions are negative for some range 
of angles, which does not happen for the temperature. To interpret the cross cor- 
relation we can consider the polarization pattern around a hot spot (T > 0). The 
cross correlation starts positive, implying a radial pattern of polarization. Not all 
the polarization is correlated with the temperature so it is hard to see this trend in 



Figure |2-1| . In Figure |2-10| we only plot the correlated part of the polarization, here 
it is clear that the vectors are preferentially radial around hot spots. 

As we move out to larger angles the cross correlation changes sign. When it 
is negative the pattern becomes tangential. For large separations the polarization 
around a hot (cold) spot is tangential (radial). A point worth noting is that the cross 
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correlation goes to zero as 9 goes to zero, in contrast to what happens for the other 
correlation functions. Symmetry arguments dictate that the Q Stokes parameter at a 
given point cannot be correlated with the temperature at that same point. What sign 
would this correlation have? Equivalently, in what direction would the polarization 
be? Only when we consider two points separated by some distance is the symmetry 
broken. The vector joining the two points becomes the privileged direction and the 
polarization can be preferentially parallel or perpendicular to this direction. 
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Chapter 3 

The Line of Sight Integration 



The field of cosmic microwave background (CMB) anisotropies lias seen a rapid devel- 
opment since its first detection by the COBE satellite only a few years ago. There are 
now several reported experimental results that are detecting anisotropies on degree 
angular scales (see 0] and for a recent review), which together with a few upper 
limits on smaller angular scales already give interesting limits on cosmological models. 
With the development of the new generation of experiments now being proposed one 
hopes to accurately map the CMB sky from arcminute scales to several degree scales. 
The amount of data thus provided would allow for an unprecedented accuracy in the 
determination of cosmological parameters. Theoretical modeling shows that CMB 
anisotropies are sensitive to most of the cosmological parameters and have a distinc- 
tive advantage over other cosmological observations in that they probe the universe 
in the linear regime. This avoids the complications caused by physical processes in 
the nonlinear regime and allows to use powerful statistical techniques to search over 
the parameter space for the best cosmological model (see e.g. 0, |, |§]). 

A large stumbling block in this program has been the speed of theoretical model 
calculations, which are still too slow to allow for a rapid search over the parameter 
space. The development of a fast and accurate algorithm to calculate the anisotropies 
becomes essential to analyze the high quality data that the future promises to deliver. 

1 Based on U. Seljak & M. Zaldarriaga, Astrophys. J. 469, 437 (1996) and M. Zaldarriaga, U. 
Seljak & E. Bertschinger, Astrophys. J. 494, 491(1998). 
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The parameter space describing the possible models is so large that a search over 
parameter space would be impossible without a fast algorithm. Very fast approximate 
methods have been developed but can only provide a 10% accuracy. In this chapter we 
present an algorithm that is both very fast (more than two orders of magnitude faster 
that the usual Boltzmann approach) and exact within linear perturbation theory. 



3.1 Einstein and Fluid Equations 

In this section we present the Einstein and fluid differential equations for the metric, 
cold dark matter (CDM) and baryons that must be solved to calculate the CMB 
anisotropy spectra produced by density perturbations in our Universe. We will also 
present the Einstein equations for gravity waves. These equations are the basis of the 
traditional methods and are also used in the integral method, discussed this chapter. 
The derivation of the Einstein and fluid equations can be found in the literature (e.g. 
0), so we only present the final results. We restrict the treatment to spatially flat 
universes. The interested reader can find the generalization to arbitrary Robertson- 
Walker background in ||. 
The metric is written as 

ds^ = —dt^ + a^{Sij + hij)dx^dx^ 

= a^[-dT^ + {5^j + hij)dx'dx^], (3.1) 

where a is the expansion factor, Xj the comoving coordinates and r = J dt/a the 
conformal time. We are using units in which c = 1. The space part of the unperturbed 
metric is a Kronecker delta, 6ij and hij is the metric perturbation in synchronous 
gauge 0. The metric perturbations have contributions from scalar (density) and 
tensor (gravity waves) modes, hij = hf^ + hj^. Although all observable quantities are 
identical in different gauges the computational efficiency to obtain them within a given 
accuracy is not. This criterion lead us to work in synchronous gauge. In comparison 
to the longitudinal gauge |]rU| it is about 20% more efficient with isentropic initial 
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conditions and even more so with isocurvature initial conditions, which are difficult 
to set up in the longitudinal gauge. 

Wc start by considering perturbations produced by density modes. When working 
with linear theory in a flat universe it is convenient to use Fourier modes because they 
evolve independently. These modes are the eigenfunctions of the Laplacian operator 
that we shall call G{k, x), 

V'^G{k,x) = -k'^G{k,x). (3.2) 

We expand all the perturbations in terms of G and its spatial covariant derivatives. 
For example, the metric perturbations for a single mode are given by 

hfj = ^S,jG -(h + 6v){k-'G\ij + ^5,jG), (3.3) 

where h and {h + 677) are the trace and traceless part of the metric perturbation. 
The perturbed Einstein's equations result in the following equations for h and rj 
(Bertschinger 1996), 

k^r]---h = -SnGa^Sp 
2 a 

k^T] = 4:nGa^{p + p)kv. (3.4) 

{G here stands for the gravitational constant and should not be confused with the 
mode functions G{k, x)); 6p and v characterize the density and velocity perturbations 
{v = ik ■ v), 5p = J2j PjSj, {p + P)v = J2j{Pj + Pj)'Vj, where pj and pj are the mean 
density and pressure of the j-th species and the sum is carried out over all the different 
species in the universe. 

The equation for the cold dark matter density perturbation 5c is, 

Sc = (3.5) 
where by definition in this gauge the cold dark matter particles have zero peculiar ve- 
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locities. The Euler equation for the baryons has additional terms caused by Thomson 
scattering and pressure. Baryons have velocities relative to the dark matter, 

6b = -kvb - ^ , 

Vb = --Vb + clk6b + ^^aneXeUriv^ - Vb) ■ (3.6) 
a 3pb 

Here Cs is the baryon sound speed, Vb is the baryon velocity, v-y is given by the 
temperature dipole v^ = 3 An and p^, pb are the mean photon and baryon densities 
respectively. The Thomson scattering cross section is 0"^, Ue is the electron density 
and Xe is the ionization fraction. 

There are two independent degrees of freedom or polarizations for gravity waves. 
The purturbed metric for each is, 

hi> ' = {ei± ie2)i (ei ± ie2)j ht{T) G{k, x), (3.7) 

the ± labels the polarization mode of the gravity wave and (ei,e2,fc) form an or- 
thogonal basis. The Einstein equations lead to 

ht + 2-ht + k'^ht = 0. (3.8) 
a 

We have neglected the source terms on the right hand side produced by the neutrino 
and photon anisotropic stress. It is straightforward to add this source, but it has 
a negligible effect in practice because for modes outside the horizon free streaming 
cannot create a shear while for modes inside the horizon the k'^ht term in (|3.8|) becomes 
dominant over the source terms. 

3.2 Boltzmann equation 

The photons are described using their distribution function which depends on time, 
spatial position and direction of propagation of the photons (n). In fact polarization 
depends also on the two axes perpendicular to n used to define the Stokes parameters. 
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The evolution of the distribution function is described using the Boltzmann equation 
which states that when one follows a light ray the temperature or polarization of the 
radiation can change for two independent reasons: gravitational redshifts or blueshifts 
and Thomson scatterings. We can then write, 

^ Qx' -^gravity -^thomsoni l"^'^/ 

X = Tor{Q ± ill) which we introduced in Chapter ^ to characterize the fractional 
fluctuations in the CMB field (ie. T = ^ = The temperature of the CMB 
blackbody can be used instead of the intensity at a particular wavelength because 
both the gravitational redshift and Thomson scattering preserve the blackbody nature 
of the spectrum. The gravity source only acts on the temperature equation and not 
on the polarization, 

Tgravity — ^ nPhij. (3.10) 

This equation is valid for both density modes and gravity waves. We will derive the 
Thomson scattering source in p.2.1| . 

3.2.1 Thomson scattering 

The Thomson scattering cross section is, 

da 3(Jt 2 



dVL Stt 



~e-~e'\\ (3.11) 



where ctt is the Thomson scattering cross section and e and e' are the unit vectors that 
describe the polarization of the electric field of the scattered and incoming radiation 
respectively. The scattering terms in equations ( |3.9| ) are most easily computed in 
the coordinate system where the incident photons travel along the z axes and the 
electrons are at rest. If fi is the direction of the incident photon and n that of the 
scattered one then n' = i = (^ = 0, = 0) and (6', cf)) describe n. For a given 
scattering event, the Thomson scattering matrix is the simplest when expressed in 
terms of the intensities of radiation parallel (Ty) and perpendicular {T±) to the plane 
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containing both n and fi . Equation 
incoming and scattered radiation, 



( p. 11 ) leads to the following relation between 



T|| 
U 



IGtt 
3 

IQn 
3 

IGtt 



cos^ 9 % 



cos 9 U'. 



(3.12) 



We have normalized the equation ( |3.12| ) so that the number of photons is conserved 
by the scattering. In this way we only need to introduce the scattering rate to obtain 
the terms needed in the Boltzmann eqaution. The total intensity is the sum of the 
two components, T = T\\ + T±, while the difference gives polarization Q = T\\ — T±. 
Because the components are measured using this coordinate system, defined by the 
plane that contains both h and n, the Stokes parameters of the incoming radiation 
Q' and U' depend on the angle (p of the scattered photon, while Q and U are already 
measured relative to the correct frame. It is more useful to refer the Stokes parameters 
of the incoming radiation relative to a fixed frame. To achieve this we construct the 
scattering matrix in terms of T', Q' + iU' = exp{2i(f)){Q' + iU') and Q' — iU' = 
exp{—2i(j)){Q' — ilf'), where we have used the transformation law (equation |2.1| ) to 
relate the two sets of Stokes parameters. 

Equation ( |3.12| ) implies that the scattered radiation in direction n which initially 
came in direction n is 



6T{h',h) = — 



1 

An 
3 



^(1 + cos' 9)T' + |(cos2 9 - l)e'^''^{Q' + tU') + 



l{cos' 9 - l)e''%Q' - lU' 



6{Q±iU){h\n) = — 



1 
3 



^(cos' 9 - l)T' + l{l± cos9fe-^"^{Q' + iU')+ 
4 8 



lTcos9fe^'^{Q'-iU'^ 



(3.13) 



We introduced a 5 to indicate that the final expression for the scattered field is an 
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integral over all directions n', 



Xifi) 



Thomson 



-aaTTle Xe 



X{n) + / dn'SX{n,h) 



(3.14) 



X stands for T and {Q ± iU). We introduced the scattering rate aarUe Xe where 
TT-e Xe is the density of ionized electrons. The first term accounts for the photons that 
are scattered away from the line of sight and the expansion factor a is introduced 
because we are calculating the derivative with respect to conformal time. 

Equation ( |3.13| ) for the scattering matrix is written in the frame where n' = {9' = 
0,0' = 0). The first spherical harmonics are explicitly ±2^2^(''^) = \/Y|^(l-Fcos^)^e^*''^, 
±2Y2'^ih) = ^3^(1 ± cosOye-^''^ and ±2Y^ih) = ^^(1 - cos^ 9) . These together 
with oYQ^^in) = 



^5^0, oYrin) = ^J^6mo, and ±2Y^{n') = ^6^^2 enable us 
to rewrite ( p.l3| ) in a more useful form {6ij is the Kronecker delta), 

6T{h',h) = aTY.[{^oYr{n)oYrin)+ oYaf^)oY,-{n'))T' 
-^/|oFr(n) 2Yr{n) {Q' + ^U') 



20 V 3 



,Yr{n) .2Yr{n'){Q' -lU') 



5{Q±iU){n\h) = (TtE[ 



6 2 
20 V 3 



^2Yr{n) ,Y^{n)r 



6 



+ 20 2Yr{n'){Q' + iU') + 

^±2rr(n) -2Y^{h'){Q'-iU'] 



(3.15) 



This form has the advantage of being independent of the coordinate system. We 
will use the scattering matrix in the frame where || i; and not n = z. Here k is 
the wavevector of the Fourier mode under consideration. The addition theorem for 
the spin harmonics ( p.l5|) can be used to show that the sum J2m sYi"^{n) s'Y"''{ri') 
acquires a phase change under rotation of the coordinate system that exactly cancels 
the phase change in the transformation of {Q ± iU) in equation ( p. 151) . We may 
therefore use this equation in any coordinate system. 
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3.2.2 Density fluctuations 

As we are dealing with a linear problem with a spatially constant background we 
may consider only one eigenmode of the Laplacian at a time. We may choose without 
loss of generahty that fe || i;. To define the Stokes parameters we use the spherical 
coordinate unit vectors (e^, e^). 

The density field produced by a single mode has two important symmetries. It is 
invariant under rotations around k and parity operations where the x or y axes change 
sign (fe II z). By invariant we mean that the transformed density field p'{x) — p{x), 
there is no prime in the argument of p'. What is directly relevant to the calculation 
of the anisotropics is that the gravity source term satisfies these symme- 

tries. The rotational symmetry implies that the neither the temperature nor the 
Stokes parameters (in the spherical basis) can depend on 0. The parity symme- 
try has interesting consequences on the allowed polarization patterns as well. For 
simplicity we focus on the Stokes parameters at the origin. Symmetry implies that 
{Q ± iU)'{h) = {Q ± iU){h) for both rotation and parity operations. We consider 
a parity operation that changes the sign of the y axes. If we consider the direction 
n = {9, 0) when we apply the transformation its location remains unchanged, h' — h. 
On the other hand e<^, used to define the Stokes parameters changes sign. As a con- 
sequence U changes sign, U'{h') — —U{h). This is inconsistent with the symmetry 
statement unless U is zero. If U is zero in this particular direction it has to be zero in 
all directions because of the rotational symmetry. In this particular coordinate sys- 
tem only Q is different from zero and we denote it by Ap, so that Ap — Q — Q ±iU 
{U — 0). The temperature anisotropy for the single eigenmode is denoted by Af. 
We introduced a A in the notation to enphasize that these are the contributions of 
one mode. The total field T, Q and U are obtained by adding the contribution of all 
modes. 

For a plane wave, rotational symmetry implies that both Af. and Ap depend only 
on the angle between n and z {k \\ z), so only harmonics with m — (which do not 
depend on 0), are needed in the expansion. To calculate the evolution of these two 
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quantities we expand them as, 



I 
I 

G{k,x)J2{ 

I 

G{k,x)J2{ 

I 

G{k,x)J2{ 



— z 



47r(2/ + 1)A|, Y,%n) 
(2/ + l)Af,PK/i) 

^47r(2/ + l)(/ + 2)!/(/-2)! sA^, 2Yi\h) 



47r(2/ + + 2)!/(/ - 2)! .sAf^ -2l^°(n) 
(2/ + 1) ±2A|,Pf (/x) (3.16) 



where G{k, x) = exp{ik ■ x) and ft = k ■ n. We added a subindex ±2 to ±2Ap; to 
denote that they are the expansion coefficients in spin ±2 harmonics [| and we used 
the exphcit expression for spin s harmonics with m = to write them in terms of 
associated Legendre polynomials ||Tl[| , 



'(2/ + i; 

Air 



Pi{cos9) 



\\ 47r (/ + 2)! ' ^ 



(3.17) 



As stated above, scalar modes in this reference frame have U = 0, so Ap describes 
both spin ±2 quantities. For m = one has = -2^° and so 2Ap; = _2Ap;. This 
is a very important result, because the fact that both coefficients are equal implies 
that this single Fourier mode does not produce any B mode polarization (equation 
|2.6|) . The observed polarization field is the superposition of that produced by each 
Fourier mode, and thus we have shown that the pattern of polarization produced by 
density perturbations has no B component. Note that B is invariant under rotations 
so the fact that we have proved this result in a particular reference frame that depends 
on k is not important. This conclusion does not apply to U: although it is zero in 
this particular frame it is not in others. When we superimpose the perturbations of 



^The relation between these coefficient and those used in Zaldarriaga & Seljak (1997) is ±2^pi = 



Va-2)!/(/ + 2)!A| 



El- 
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all modes B remains zero while U does not. 

We can replace equation ( |3.1(j| ) into ( |3.9| ) and use the expresion for the gravita- 
tional redshift, equations (|3.3|) and ( |3.10|) , to obtain the Boltzmann equation for the 
CMB photons; 



- ^(^ + 6?7)^2(/i) + ^TlThomson 



1 

3' 



^P\Thomson- 



(3.18) 



The first term in the temperature equation represents the effect of gravitational red- 
shift, while A^T\Thomson and p\Thomson are the changes in the photon distribution 
function produced by Thomson scattering. After inserting equation ( |3.16|) into equa- 
tion ( p. 18 ) one obtains a system of two coupled hierarchies, one for the temperature 
and the other for polarization, 



A^ 



A^ 



A^ 



^ + A^ 

6 



TO\Thomson 



+ A^ , 



2 

15 



+ —k^a + A'^ 



T2\Thomson 



k 



21 + 1 
k 

21 + 1 



l^T(i-i) - + l)^r(/+i)] + A. 



(/ - 2) 2A^,_, - (/ + 3) 2A^,^i 



Tl\Thomson 

+ sA^ 



J > 2 



2'-^Pl\Thomson'> 



(3.19) 



where a = {h + 6fi)/2k'^, and we wrote sepatarely the first 3 equations of the tempera- 
ture hierarchy that contain the gravity sources. We also used the recurrence relations 
for the Legendre functions. 



21 + 1 



[/Pz_i + (/ + l)P, 



21+ 1 



{I + 2)Pl, + {l- l)Pl, 



(3.20) 



Substituting the expansion for the Stokes parameters from equation (|3.16|) into 
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equation (|3.15| ) and using equation (|3.14|) we find 



Tl I Thomson 



I 

±2'-^Pl\ Thomson 



A^i + J dQoYi"'{h)ST{n) 



n 



-aaTXeUe 



±2Af, + I dn^Yl^ifi) 6{Q±tU)ih) 



n 

20' 



(3.21) 



with 



n = A 



T2 



6( 2A 



s 

P2 



^2A' 



P2) 



(3.22) 



The differential optical depth for Thomson scattering is denoted k = an^XeCrT- Note 
that the polarization has sources only at / = 2. Equation ( p.21| ) is valid in the rest 
frame of the electrons, so in the reference frame where the baryon velocity is the 
distribution of scattered radiation has an additional dipole. The final expression for 
the Boltzmann hierarchy is 
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The last important ingredient is the relation between A^ p-^ for each Fourier mode 
and the observed CMB power spectra. For a single mode we can combine equations 
(g^D and (|316| ) to show 
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with A|_, = -^(/ + 2)!/(/-2)! ±2A|,. 

The different Fourier modes are statistically independent so the observed total 
power in the integral of the powers produced by all modes, 

Ci^T,E)l = {A7ryJedkP{k)\Afr,E)l{^,T = To)\' 

Cci = (47r)2 1 edkP{k)A-^i{k, T = ro)A|;(A;, r = tq) (3.25) 

with P{k) the primordial power spectrum. The primordial power spectrum gives the 
power in each mode at the initial time. For example if we choose to normalize the 
initial conditions so that 6c = 1 then P{k) gives the initial power in the cold dark 
matter component. 

3.2.3 Gravity waves 

For the gravity wave modes the source term in the temperature equation is | fi^fi^hj^^^^ - 
I (1 — cos^ 6')e^*^'^ ht G{k,x). There is an explicit (j) dependence of the source so 
the symmetries of rotation and parity satisfied by the density modes discussed in 
the previous section no longer apply. This implies that U is no longer zero and that 
gravity waves will produce a non-zero B. This is a crucial point and may allow in 
the future to detect the stochastic background of gravity waves if it is present in our 
universe. We will discuss this further in section |4.3| . 
For gravity waves we can expand, 
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X _2A2^) -2Y,^\fi). (3.26) 
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Because of the dependence of the source term we had to use harmonics with m — ±2 
for the two different gravity wave polarizations. 

Now that PI 7^ -2^ PI , we need two different hierarchies to solve for both 
of them separately. The derivation of the hierarchies is analogous to that for the 
density modes yielding, 
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We can change variables to obtain hierarchies for E and B directly, 
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with n = AJ^^^ - 6(+2AJ^^^ + -sAJ^^^) = (Ajf ^ + V6A^f We introduced 

tildes in the notation to remind the reader that we introduced extra I dependent fac- 
tors for later convenience. Free streaming couples E and B modes, but the Thomson 
scattering only creates E (the kIl/20 term is only present in the E hierarchy). It is 
through freestreaming and not scattering that B is created. 

Finally the relation between the observed CMB power spectra and the perturba- 
tions for each mode is, 

C^T,E,B)i = i4n)'lk'dkPtik)\A^I%^^^{k,r^ro)\' 

Cci = {inf J k^dkPt{k)A'f^{k,T = To)Alf\k,T = To), (3.32) 

where Ph{k) denotes the initial power spectrum for the gravity waves and we normal- 
ized the initial conditions to ht — 1 

The results for scalars (m = 0) and tensors (m = 2) can all be combined (the 
following formulas also applies to vector modes with m — 1). We expand all variables 
as. 
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The expansion coefficients will satisfy the hierarchy, 
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with n = A^2^^ — 6(+2Ap2'^'' + _2A^^^''). Similar expressions can be written for the 
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E and B hierarchies. Equation ( 3.34 ) does not include the gravitational redshift term 
or the Doppler shift caused by the the baryon velocity which must be must be added 
to the temperature equation. In the general case we can define 
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and use equation (|3.32|) to calculate the observer power spectra. 



3.3 The Line of Sight Method 

Instead of solving the coupled system of differential equations ( |3.23| ) one may formally 
integrate equations (|3.18|) along the photon past light cone to obtain (e.g. [p!4|), 

Af(fc,ro) = / dre'^^("-"«)e-'^{/ie-''[Afo + ^/it^b + -P2(/u)n] 
Jo 2 

Al,(k,To) = ] r dre'^^^^-^'^^e-'^kllP^U). (3.36) 
4 JO 

These lead to equations for the multipoles, 

Af;(A;,r = ro) = / S^{k,T)3i{x)dT 
Jo 

,Af,i{k,T = To) = r S'p{k,T)^-^dT 

Jo 

X = k{ro-T) (3.37) 

where the source functions are, 
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S'piKT) = -^gU. (3.38) 

where a = {h + 6fi)/2k'^, ji{x) are the spherical Bessel functions and g is the visibihty 
function g = kexp{—K,). The source function for the temperature looks rather com- 
plicated because we have done several integrations by parts to express everything in 
terms of ji{x) rather than its derivatives. The polarization solution is straightforward 
to check. It is enough to substitute equation ( p.38|) into the hierarchy (|3.23| ) using 
the limit lim^.^o Ji(^)/^^ = 1/15 and the recursion, 



d Jijx) ^ 1 
dx^ x2 ' (2/ + 1) 



(3.39) 



Some of the terms in the source function Sj<{t) are easily recognizable. It may be 
simpler to express them in the more familiar conformal Newtonian gauge, Afp^'^^+a = 
A'^Q^ + vl^^*^^ + ka = vl°^^ and ci + ?) = </) + ■?/'. Here (j) and ijj are the conformal 
gauge gravitational potentials. All moments higher than / = 1 are gauge invariant, so 
in fact the source terms S(t,p) must be gauge invariant. The first two contributions in 
the source Sj- (equation p. 38]) A^o + a are the intrinsic temperature and gravitational 
potential contributions from the last-scattering surface, while the third contribution 
is part of the velocity term, the other part being k~^gvb + a in the next line. These 
terms make a dominant contribution to the anisotropy in the standard recombination 
models. The first term in the second line is the so-called integrated Sachs- Wolfe 
term and is important after recombination. It is especially important if matter- 
radiation equality occurs close to the recombination or in i^mattcr 7^ 1 models. In both 
cases the gravitational potential decays with time, which leads to an enhancement 
of anisotropics on large angular scales. Finally we have the terms caused by photon 
polarization and anisotropic stress, which contribute to 11. These terms affect the 
anisotropy spectra at the 10% level and are important for accurate model predictions. 
Moreover, they are the sources for photon polarization. 

The main advantage of (|3.38|) is that it decomposes the anisotropy into a source 
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term p, which does not depend on the multipole moment / and a geometrical 
term ji, which does not depend on the particular cosmological model. The latter only 
needs to be computed once and can be stored for subsequent calculations. The source 
term is the same for all multipole moments and only depends on a small number of 
contributors in ( ^.38| ) (gravitational potentials, baryon velocity and photon moments 
up to / = 2). By specifying the source term as a function of time one can compute 
the corresponding spectrum of anisotropics. Equation ( p.38| ) is formally an integral 
system of equations, because the first moments appear on both sides of equations. 
To solve for these moments it is best to use the equations in their differential form 
( p.23| ), instead of the integral form above. Once the moments that enter into the 
source function are computed, one can solve for the higher moments by performing 
the integration in (|3.38|) . 

The solution for the tensor modes can similarly be written as an integral over the 
source term. 
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Equations ( p.38| ) and ( p.4(]| ) are the main equations of this chapter and form the 
basis of the line of sight integration method of computing CMB anisotropics. The 
treatment can be generalized to arbitrary FRW backgrounds 0. In the next section 
we will discuss in more detail the computational advantages of this formulation of the 
Boltzmann equation and its implementation. 



64 



3.4 Calculational Techniques 



In the previous section we presented the expressions needed for the implementation 
of the hne of sight integration method. As shown in equations (|3.38| ) and (|3.40|) one 
needs to integrate over time the source term at time r multiphed with the sperical 
Bessel function evaluated at x = k{TQ — r). The latter does not depend on the 
model and can be precomputed in advance. Fast algorithms exist which can compute 



spherical Bessel functions on a grid in x and I in short amount of time (e.g. ||16|| ). 
The grid is then stored on a disk and used for all the subsequent calculations. This 
leaves us with the task of accurately calculating the source term, which determines 
the CMB spectrum for a given model. Below we discuss some of the calculational 
techniques needed for the implementation of the method. We especially highlight the 
differences between this approach and the standard Boltzmann integration approach. 
Our goal is to develop a method which is accurate to 1% in Ci up to / ~ 1000 over 
the whole range of cosmo logical parameters of interest. These include models with 
varying amount of dark matter, baryonic matter, Hubble constant, vacuum energy, 
neutrino mass, shape of initial spectrum of perturbations, reionization and tensor 
modes. The choice of accuracy is based on estimates of observational accuracies that 
will be achievable in the next generation of experiments and also on the theoretical 
limitations of model predictions (e.g. cosmic variance, second order effects, etc.). 
Most of the figures where we discuss the choice of parameters are calculated for 
the standard CDM model. This model is a reasonable choice in the sense that it 
is a model which exhibits most of the physical effects in realistic models, including 
acoustic oscillations, early-time integrated Sachs- Wolfe effect and Silk damping. One 
has to be careful however not to tune the parameters based on a single model. We 
compared our results with results from other groups [p!7| , p!8[] for a number of different 
models. We find a better than 1% agreement with these calculations over most of the 
parameter space of models. The computational parameters we recommend below are 
based on this more detailed comparison and are typically more stringent than what 
one would find based on the comparison with the standard CDM model only. 
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Figure 3-1: CMB spectra produced by varying the number of evolved photon mul- 
tipole moments, together with the relative error (in %) compared to the exact case. 
While using = 5 produces up to 2% error, using = 7 gives results almost identical 
to the exact case. 

3.4.1 Number of coupled differential equations 

In the standard Boltzmann method the photon distribution function is expanded to 
a high /max ( ^.23| ) and typically one has to solve a coupled system of several thousand 
differential equations. In the integral method one evaluates the source terms S{k, r) 
as a function of time (equations p.38| and |3.4CI| ) and one only requires the knowledge 



of photon multipole moments up to / = 2, plus the metric perturbations and baryon 
velocity. This greatly reduces the number of coupled differential equations that are 
needed to be solved. For an accurate evaluation of the lowest multipoles in the integral 
method one has to extend the hierarchy somewhat beyond 1 = 2, because the lower 
multipole moments are coupled to the higher multipoles ( p.23| ). Because power is only 
being transferred from lower to higher / it suffices to keep a few moments to achieve a 
high numerical accuracy of the first few moments. One has to be careful however to 
avoid unwanted reflections of the power being transferred from low / to high /, which 
occur for example if a simple cut-off in the hierarchy is imposed. This can be achieved 
by modifying the boundary condition for the last term in the hierarchy using the free 
streaming approximation |12|, [T^. In the absence of scattering (the so-called free 



streaming regime), the recurrence relation among the photon multipoles in equation 
( ^.23| ) becomes the generator of spherical Bessel functions. One can therefore use 
a different recurrence relation among the spherical Bessel functions to approximate 
the last term in the hierarchy without reference to the higher terms. The same 
approximation can also be used for polarization and neutrino hierarchies. This type 
of closure scheme works extremely well and only a few multipoles beyond I = 2 are 



needed for an accurate calculation of the source term. This is shown in Figure 3-1 



where a relative error in the spectrum is plotted for several choices of maximal number 
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Figure 3-2: Relative error between the exact and interpolated spectrum, where 
every 20th, 50th or 70th multipole is calculated. The maximal error for the three 
approximations is less than 0.2, 0.4 and 1.2%, respectively. The rms deviation from 
the exact spectrum is further improved by finer sampling, because the interpolated 
spectra are exact in the sampled points. For the sampling in every 50th multipole 
the rms error is 0.1%. 

of photon multipoles. We choose to end the photon hierarchy (both anisotropy and 
polarization) at = 8 and massless neutrino at 1^ = 7, which results in an error 
lower than 0.1% compared to the exact case. Instead of a few thousand coupled 
differential equations we therefore evolve about 35 equations and the integration time 
is correspondingly reduced. 



3.4.2 Sampling of CMB multipoles 

In the standard Boltzmann integration method one solves for the whole photon hi- 
erarchy ( 3.23| ) and the resultant A; is automatically obtained for each / up to some 



^max- The CMB spectra are however very smooth (see Figure |3-1|) , except for the 
lowest / where the discrete nature of the spectrum becomes important. This means 
that the spectrum need not be sampled for each / and instead it suffices to sparsely 
sample the spectrum in a number of points and interpolate between them. Figure |3-2| 
shows the result of such interpolation with cubic splines (see e.g. [|T^) when every 
20th, 50th or 70th / is sampled beyond / = 100 with an increasingly denser sampling 
towards small /, so that each / is sampled below / = 10. While sampling of every 70th 
/ results in maximal error of 1%, sampling in every 20th or 50th / gives errors below 
0.2 and 0.4%, respectively. We choose to compute every 50th Ci beyond I = 100 in 
addition to 15 / modes below I = 100, so that a total of 45 / modes are calculated 
up to Zmax = 1500. This gives a typical (rms) error of 0.1%, with excursions of up to 
0.4%. The number of integrals in equation ( ^.25] ) is thus reduced by a factor of 50 at 
high / and the computational time needed for the integrals becomes comparable to 
or smaller than the time needed to solve for the system of differential equations. 
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Figure 3-3: Error in the spectrum caused by insufficient temporal sampling of the 
source term. Inaccurate sampling of the source during recombination leads to numer- 
ical errors, which can reach the level of 1% if the source is sampled in only 10 points 
across the recombination epoch. Finer sampling in time gives much smaller errors 
for this model. Comparisons with other models indicate that sampling in 40 points 
is needed for accurate integration. 

3.4.3 Integration over time 



For each Fourier mode k the source term is integrated over time r (equation 3.38| ). The 



sampling in time need not be uniform, because the dominant contribution arises from 
the epoch of recombination around z ~ 1100, the width of which is determined by the 
visibility function g and is rather narrow in look-back time for standard recombination 
scenarios. During this epoch the sources acoustically oscillate on a time scale of Csk~^, 
so that the longest wavelength modes are the slowest to vary. For short wavelengths 
the rate of sampling should therefore be higher. Even for long wavelengths the source 
function should still be sampled in several points across the last-scattering surface. 
This is because the terms in ( ^.38 ) depend on the derivatives of the visibility function. 



If the visibility function g is narrow then its derivative will also be narrow and will 
sharply change sign at the peak of g. Its integration will lead to numerical roundoff 
errors if not properly sampled, even though positive and negative contributions nearly 
cancel out when integrated over time and make only a small contribution to the 



integral. Figure 3^ shows the error in integration caused by sampling this epoch 
with 10, 20 or 40 points. Based on comparison with several models we choose to 
sample the recombination epoch with 40 points, which results in very small (~ 0.1%) 
errors. After this epoch the main contribution to the anisotropics arises from the 
integrated Sachs- Wolfe term. This is typically a slowly changing function and it is 
sufficient to sample the entire range in time until the present in 40 points. The 
exceptions here are models with reionization, where the visibility function becomes 
non-negligible again and a new last-scattering surface is created. In this case a more 
accurate sampling of the source is also needed at lower redshifts. 
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Figure 3-4: In (a) Af^^5Q(A;) is plotted as a function of wavevector k. In (b) A|^ ]^5q(A;) 
is decomposed into the source term Sj^ integrated over time and the spherical Bessel 
function ji5o(fcTo)- The high frequency oscillations of j^5q(/c) are caused by oscilla- 
tions of the spherical Bessel function ji5o(^To), whereas the source term varies much 
more slowly. This allows one to reduce the number of k evaluations in the line of 
sight integration method, because only the source term needs to be sampled. 

3.4.4 Integration over wavenumbers 

The main computational cost of standard CMB calculations is solving the coupled sys- 
tem of differential equations. The number of /c-modes for which the system is solved is 
the main factor that determines the speed of the method. For results accurate to /max 
one has to sample the wavenumbers up to a maximum value /cmax = ^max/To- In the 
line of sight integration method solving the coupled system of differential equations 
still dominates the computational time (although for each mode the time is signifi- 
cantly shorter than in the standard Boltzmann method because of a smaller system of 
equations). It is therefore instructive to compare the number of k evaluations needed 
in each of the methods to achieve a given accuracy in the final spectrum. 

In the standard Boltzmann method one solves for A'^i{k) directly, so this quantity 
must be sampled densely enough for accurate integration. Figure |3-4| a shows Af ^(A;) 
for / = 150 in a standard CDM model. One can see that it is a rapidly oscillating 
function with a frequency k ~ Tq^. Each oscillation needs to be sampled in at least a 
few points to assure an accurate integration. To obtain a smooth CMB spectrum one 
typically requires 6 points over one period, implying 2/max k-mode evaluations. This 
number can be reduced somewhat by filtering out the sampling noise in the spectrum 
p!9| , but even in this case one requires at least 1-2 points per each period or /max/2 
fc-mode evaluations. 

To understand the nature of these rapid oscillations in Af ^(A;) we will consider 
wavelengths larger than the width of the last scattering surface. In this case the Bessel 
function in ( |3.38D can be pulled out of the integral as jiikro) because the time at which 



recombination occurs, when the dominant contribution to A|.;(/c) is created, is much 
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smaller than tq and kAr ^ 1 (Ar is the interval of time for which the visibility 
function differs appreciably from zero). So the final ^(/c) is approximately the 
product of ji{kTo) and 5*1^ integrated over time, if the finite width of the last scattering 
surface and contributions after recombination can be ignored. 

Figure ^-4] b shows the source term 5*^ integrated over time and the Bessel function 
ji{kTQ). It shows that the high frequency oscillations in A^i{k) seen in Figure P-4| a 
are caused by the oscillation of the spherical Bessel functions, while the oscillations 
of the source term have a much longer period in k. The different periods of the 
two oscillations can be understood using the tight coupling approximation |jT9|, ^ 



Prior and during recombination photons are coupled to the baryons and the two 
oscillate together with a typical acoustic timescale ~ Tj-cc/ ~ tq/ VS-^rcc ~ to/50. 
The frequency of acoustic oscillations k ~ is therefore 50 times higher than the 
frequency of oscillations in spherical Bessel functions, which oscillate as Tq^. 

Because an accurate sampling of the source term requires only a few points over 
each acoustic oscillation, the total number of k evaluations in the integral method can 
be significantly reduced compared to the standard methods. Typically a few dozen 
evaluations are needed over the entire range of k, compared to about 500 evaluations 
in the standard method when a noise filtering technique is used and 2000 otherwise 
(for /max ~ 1000). Once the source term is evaluated at these points one can inter- 
polate it at points with preevaluated spherical Bessel functions, which can be much 
more densely sampled at no additional computational cost. The end result is the 
same accuracy as in the standard method, provided that the source is sampled in 
sufficient number of points. Figure p-5| shows the relative error in the CMB spectrum 
for the cases where the source term is calculated in 40, 60 and 80 points between 
and kro = 3000 (for /max = 1500). While using 40 points results in up to 1% errors, 
using 60 points decreases the maximum error to below 0.2% for this model. In general 
it suffices to use /max/30 k modes, which is at least an order of magnitude smaller 
than in the standard methods. Note that with this method there is no need to filter 
the spectrum to reduce the sampling noise, because the latter is mainly caused by 
insufficient sampling of the spherical Bessel functions, which are easy to precompute. 
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Figure 3-5: Error in the spectrum caused by insufficient k-mode sampling of the 
source term. Samphng the source with 40 points up to k = 2/inax leads to 1% errors, 
while with 60 or 80 points the maximal error decreases to 0.2%. Comparisons with 
other models indicate that sampling in 60 points is sufficient for accurate integration. 

The additional operations needed for a higher sampling (summation and source inter- 
polation) do not significantly affect the overall computational time. Moreover, if each 
Ci is accurately calculated they can be sparsely sampled and interpolated (section 
3.2). This would not be possible if they had a significant noise component added to 
them. 

3.5 Summary 

In this chapter we presented a new method for accurate calculations of CMB anisotropy 
and polarization spectra. The method is not based on any approximations and is 
an alternative to the standard Boltzmann calculations, which are based on solving 
large numbers of differential equations. The approach proposed here uses a hybrid 
integro-differential approach in solving the same system of equations. By rewriting 
the Boltzmann equations in the integral form the solution for the photon anisotropy 
spectrum can be written as an integral over a source and a geometrical term. The ffist 
is determined by a small number of contributors to the photon equations of motion 
and the second is given by the radial eigenfunctions, which do not depend on the 
particular cosmological model, but only on the geometry of space. 

One advantage of the split between geometrical and dynamical terms is that it 
clarifies their different contributions to the final spectrum. A good example of this is 
the temperature anisotropy in the non-fiat universe, which can be be written using 
a similar decomposition, except that spherical Bessel functions have to be replaced 
with their appropriate generalization [|^. This is discussed in more detail elsewhere 
0. Here we simply remark that replacing radial eigenfunctions in a non-fiat space 
with their fiat space counterpart (keeping comoving angular distance to the LSS 
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unchanged) is only approximate and does not become exact even in the large / (small 
angle) limit. The geometry of the universe leaves its signature in the CMB spectra in 
a rather nontrivial way and does not lead only to a simple rescaling of the spectrum 
by n-^L i- 

The main advantage of our line of sight integration method is its speed and ac- 
curacy. For a given set of parameters it is two orders of magnitude faster than the 
standard Boltzmann methods, while preserving the same accuracy. We compared our 
results with the results by Sugiyama and by Bode & Bertschinger |T^. In both 
cases the agreement was better than 1% up to a very high / for all of the models we 
compared to. 

The method is useful for fast and accurate normalizations of density power spec- 
tra from CMB measurements, which for a given model require the CMB anisotropy 
spectrum and matter transfer function, both of which are provided by the output of 
the method. Speed and accuracy are even more important for accurate determination 
of cosmological parameters from CMB measurements. In such applications one wants 
to perform a search over a large parameter space, which typically requires calculating 
the spectra of a several thousand models (e.g. 0). One feasible way to do so is to use 
approximation methods mentioned in the introduction. These can be made extremely 
fast, but at a cost of sacrificing the accuracy. While several percent accuracy is suf- 
ficient for analyzing the present day experiments, it will not satisfy the requirements 
for the future all-sky surveys of microwave sky. Provided that foreground contribu- 
tions can be succesfully filtered out one can hope for accuracies on the spectrum close 
to the cosmic variance limit, which for a broad band averages can indeed reach below 
1% at Z > 100. It is at this stage that fast and accurate CMB calculations such as the 
one presented in this chapter will become crucial and might enable one to determine 
many cosmological parameters with an unprecedented accuracy. 
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Chapter 4 

Information in the CMB 
Polarization^ 

4.1 Introduction 

It is now well established that temperature anisotropies in CMB offer one of the best 
probes of early universe, which could potentially lead to a precise determination of 
a large number of cosmological parameters |jl|, ^, ^. The main advantage of CMB 
versus more local probes of large-scale structure is that the fluctuations were created 
at an epoch when the universe was still in a linear regime. While this fact has long 
been emphasized for temperature anisotropies, the same holds also for polarization in 
CMB and as such it offers the same advantages as the temperature anisotropies in the 
determination of cosmological parameters. The main limitation of polarization is that 
it is predicted to be small: theoretical calculations show that CMB will be polarized 
at the 5-10% level on small angular scales and much less than that on large angular 
scales @]. However, future CMB satellite missions (MAP and Planck) will be so 
sensitive that even such small signals will be measurable. Even if polarization by itself 
cannot compete with the temperature anisotropies, a combination of the two could 
result in a much more accurate determination of certain cosmological parameters, 

^Based on M. Zaldarriaga, Phys. Rev. D 55, 1822 (1997), U. Seljak & M. Zaldarriaga, Phys. 
Rev. Lett. 78, 2054 (1997) and D. N. Spergel & M. Zaldarriaga, Phys. Rev. Lett. 79, 2180 (1997). 
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in particular those that are hmited by a finite number of multipoles in the sky (i.e. 
cosmic variance). 

4.2 The Imprint of an Early Reionization of the 
Universe 

It has been pointed out that an early reionization of the universe will greatly enhance 
polarization [Q. The fact that in universes that never recombined the polarization 
would also be large was noted in many of the above studies. More recently Ng & 
Ng p discussed the polarization generated in reionized universes with instantaneous 
recombination. The Sachs- Wolfe effect was the only source of anisotropics that they 
included. They concluded that reionization at a moderate redshift could boost po- 
larization to the level of a few percent of the temperature perturbations. To make 
detailed predictions for an experiments (such as that being built at Wisconsin Uni- 
versity) a realistic recombination history should be used since polarization is very 
sensitive to the duration of recombination ||^, Baryons should also be included 
in the calculation as the acoustic oscillation in the photon-baryon plasma are very 
important to determine polarization. 

In this section we discuss in detail the physics behind the polarization generated 
in models where there was an early reionization after the usual recombination. These 
models show very distinct features in the polarization power spectrum including a 
new peak at low /. This peak is not present either in the standard recombination 
scenarios or in the cases where the universe never recombined and it is the cause for 
the boost in the polarization. 

4.2.1 Standard Ionization History 

In this section we review previous results for the CMB polarization for a standard 
ionization history in a flat space-time. The Boltzmann equations for the perturbations 
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Figure 4-1: /(/ + l)Ci/27r for both temperature (a) and polarization (b) for standard 
CDM and a model where the optical depth to recombination is Uri = 1-0. 



in the scalar case are given by equations (|3.18|) and ( |3.21|) , 



At + ikfiAx 
Ap + ikfiAp 



O 6 



-Ap + -p|(^)n 



-At + Ato + if^Vb + -P2(/i)n 



(4.1) 



They can be formally integrated to give (|3.3 



Arik, To) = r dre^'^'^^^-^'^h-^ike-^lATo + tfiVf, + ^P2(;u)n] 
JO 2 



1 



Ap(k,To) = - r dTe'^^^^^^'''^e-''kIlPiU). 
4 JO 



(4.2) 



Figure |4-1| shows the temperature and polarization Ci spectra for the standard 
CDM model {VLq = 1, Hq = 50kmsec~^ and Qb = 0.05), normalizing the result to 
the COBE measurement. Normalization was carried out using the fits to the shape 
and amplitude of the 4 year COBE data described in which approximately fixes 
Ctw ~ 4:7 fiK^. In Sections |2.2.2| and |2]^ we discussed an approximate solution for 
wavelengths longer than the width of the last scattering surface, Atd (equation |2.44| ). 
We showed that the polarization perturbation is 0, 



Ap = 0.51(1 -/i2)e^Mri,-ro)^^^^^^^(^^^ 



(4.3) 



where td is the conformal time of decoupling. Note that in the tight coupling regime 
A^i oc Vb- The polarization is proportional to the velocity difference between 
places separated by a distance Atd, the distance photons travel on average during 
decoupling. For the standard adiabatic initial conditions A^i and the baryon velocity 

•^This velocity is in the longitudinal gauge, w™"-'^ = v^^"^'' + ka. 
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Figure 4-2: Visibility function for standard CDM with reionization such that the 
optical depth to recombination is n^i = 1.0. 

vanish as fcr — which together with the fcAr^) factor in the previous expression 
explains the dramatic decrease of polarization with decreasing /. For large wavelengths 
the quadrupole generated in the photon distribution as photons travel between their 
last scatterings is extremely small both due to the small distance they can travel 
compared to the wavelength as well as to the small velocity differences generated by 
these small k perturbations. 



4.2.2 The Reionized Case 

In this section we consider models with early reionization and try to explain the origin 
of the new features that appear in the polarization power spectrum. For definiteveness 
we use a standard CDM model with specified optical depth to recombination Kri. For 
example if reionization occurred at a redshift of around z^i ~ 100 then K^i = 1-0, 
assuming that all hydrogen atoms are ionized up to the present epoch {xe = 1.0). 
Figure [4.2.2| shows the visibility function, g{T) = /texp(— k), for = 1.0 . The 



visibility function has a very simple interpretation: the probability that a photon 
reaching the observer last scattered between r and r + rfr is g{T)dT. The first peak in 



Figure [4.2.2| , occurring at r 120 Mpc for sCDM {h = 0.5) accounts for the photons 
that last scattered at recombination. The area under this peak, the probability that a 
photon came directly to us from this epoch, is exp(— K^j). The area under the second 
peak gives the probability that a photon scattered after reionization and is equal to 

1 — exp( — Kri). 

Figure |4-1| shows the result of numerically integrating the Boltzmann equations 
for this reionized case. On small angular scales, the polarization acoustic peaks are 
suppressed, just like in the temperature case. This is very easy to understand: only 
a fraction exp(— /t^j) of the photons reaching the observer come from recombination, 
so their contribution to the Ci power spectrum is reduced by a factor exp(— 2Kri). On 
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large scales new temperature anisotropies are created that compensate the reduction 
and leave the spectrum approximately the same. New peaks appear in the polarization 
power spectrum while the temperature anisotropy shows no new peaks. These new 
peaks in the polarization are what boost it and may make it detectable. 

Let us try to understand the origin of these peaks. Before the decoupling of 
photons and electrons, for low values of k the largest perturbation in the photon 
distribution function is the monopole, A7-0 because of the tight coupling between 
photons and electrons before recombination. Both the dipole and the quadrupole 
as well as the polarization perturbations are much smaller. But after photons and 
electrons decouple, all the temperature multipoles can grow by free streaming. Power 
is being carried from the low multipole moments to higher ones. This is a geometrical 
effect. The temperature quadrupole is growing by free streaming after recombination 
and by the time the universe reionizes there is an appreciable quadrupole that can 
generate polarization. The structure of this quadrupole explains the new features in 
the polarization power spectrum. 

The formal line of sight solution for the polarization perturbation is 

Ap = ^ r dTe''"'^^-^'>^e-^knP^(fx). (4.4) 
4 Jo 

The visibility function has two peaks one at recombination and one at reionization. 
It is convenient to separate the previous integral in two parts, 

Ap = jP^ifx) ( T" dTe''"'^^-^°^ke-''U + r dTe''"'^^-^°^ke-^n) (4.5) 

where r^j is the conformal time of the start of reionization. The first integral represents 
the polarization generated at recombination and can easily be shown to be 

A^p^ = tPK/x) r" cire'^^("-"«)fi;e-"n = e'^-A^^ (4.6) 
4 Jo 

where Ap^ is the polarization that would be measured if there was no reionization, as 
discussed in the previous section. This contribution is damped because only a fraction 
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Figure 4-3: /(/ + l)Gp/27r (a) for CDM models with varying = 0.5, 1.0, 1.5 
and (b) for models with varying cosmological constant Q\ = 0.3, 0.5, 0.7 and 
a fixed redshift of reionization Zri = 100. Reionized (k^ = 1-0) CDM mod- 
els (c) with varying = 0.3, 0.5, 0.8 and (d) with different Hubble constants 
Hq = 60, 80, 100 km sec~^Mpc~^. In all cases reionization was assumed to be 
complete (xe = 1) 

exp{—K,ri) of the photons that arrive to the observer came directly from recombination 
without scattering again after reionization. 

Let us now consider the new contribution arising from reionization. The polar- 
ization source is 11 = At2 — 6(+2Ap2 + -2^P2)- The temperature quadrupole At2 
is large and originates in the free streaming of the monopole at recombination, while 
the polarization terms do not grow after decoupling and are thus negligible to first 
approximation. Equation (|4.5| ) shows that the new polarization is basically an aver- 
age of the value of the temperature quadrupole during reionization. This accounts 
for all the new features in the polarization power spectrum of Figure |4-1| . 

To understand the origin of these new peaks let us find the amplitude of the 
temperature quadrupole at the time reionization starts r^. The monopole at recom- 
bination is approximately given by ||TU| 



(Ato + = ^^(1 + 3i?) cos(A;c,r^) - (4.7) 

where is a constant (the value of the conformal gauge gravitational potential). 



R = 3pb/4:p^\T-jj ~ SOilbh"^ and Cg = l/y3(l + R) is the photon-baryon sound speed. 
The quadrupole at r^j arising from the free streaming of this monopole is simply 
(equation |3.38| ), 

AT2(r„) = {Ato + «)(rD)j2[A;(r„ - td)] (4.8) 

where j2 is the I = 2 spherical Bessel function. 

The peaks of At2 (Tri) as a function of k will show up in the polarization power 
spectrum. The first peak of equation ( |4.8|) is approximately at the first peak of the 
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Bessel function because CsT^ ^ {r^ — t^,). The wavenumber for this first peak is 
approximately given by k{Tri — t/j) ~ 2, translating into an / value as usual according 
to I ~ k{TQ — Tri). The physical size of the wavelength is translated to an angular 
size in the sky using the distance to reionization, tq — th- In terms of multipoles, 
the first peak is at Z ~ k{TQ — T,.i) ~ 2(ro — Tn)/{Tj.i — td) 2^^. For the case 
under consideration, z ~ 100 this means / ~ 20 which agrees very well with the the 
first peak in Figure [4-3| . Only the first few peaks appear because the reionization 
scattering surface is wide and thus for smaller wavelengths the source in equation 
( [4.5|) oscillates during its width and cancels out after integration. This cancellation 
makes the new polarization small and thus hidden under the polarization generated 
at recombination. 

The major factor determining the difference in height of these new peaks for 
different models is the fraction of photons reaching the observer that last scattered 
after reionization, l — exp{—Kri). The ratio of the distances between the observer and 
reionization to that between the two scattering surfaces determines the positions of 
the peaks, and the optical depth determines their heights. 



To further illustrate these points. Figure [4-3| a shows the Cei spectrum for standard 
CDM models with varying optical depths Kri. The peaks not only vary in height but 
also in position, as the redshift of reionization increases when does and so the 
position of the peaks gets moved to smaller scales, Ipeak ~ 2^5^. Figure P3| b on 
the other hand show how these peaks vary with the cosmological constant for a fixed 
reionization redshift z^i = 100. The positions hardly change as both the distance to 
reionization and the distance between the two scattering surfaces scale approximately 
in the same way with the matter density (in this calculations the matter density is 
given by Qq = 1—Qa where is the energy density due to the cosmological constant). 
The distance to a fixed redshift increases with the cosmological constant, the optical 
depth increases, and consequently the peaks should get higher. The fact that this 
is not the case is a consequence of the COBE normalization, because models with 
larger values of the cosmological constant have larger additional contributions to the 
low / temperature anisotropies from the ISW effect while polarization is not affected 
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by the late time variation of the gravitational potential which does not affect the 
polarization. Thus the changes in the normalization to keep the value of Ctw fixed 
partially compensate the change in the height of the new polarization peaks produced 
by the larger optical depth. 

Figures 3c and 3d explore the dependence of the polarization power spectrum with 
the baryon density and the Hubble constant for a fixed optical depth to decoupling, 
Kri = 1-0, keeping the rest of the parameters the same as sCDM. The height of the 
first peak in the spectrum remains nearly constant as it is determined by Kri which 
was kept fixed. The fact that the peaks move is due to the fact that the redshift of 
reionization is changing, it is given by (1 + Zri) ~ 100[Kr.i(0.5//;.)(0.05/f2{,)(l/a;e)]^''^ 
and so / scales approximately as / oc (Kri/hilbXeY^^- 

In the sCDM model the reionization must have occurred extremely early (z^j ~ 
100) in order to produce an optical depth of unity; even an optical depth of = 0.5 
is only obtained for Zri ~ 60. The situation is different for open models or models 
with a cosmological constant. An approximate scaling for the optical depth valid 
for ^lo^ri ^ 1 (/ififeXe/fiy^)(l + Zri)^^'^, SO for example reionization starting 

at Zri ~ 23 will produce an optical depth ~ 0.5 in a model with Qo = 0.2, 
Ho = 70km sec-i Mpc'^ and = 0.1. 

The polarization increase on large scales produced by an early reionization of the 
universe can have an important impact on the accuracy with which future satellite 
missions will be able to reconstruct cosmological parameters from CMB. This will 
be the focus of chapter ^. There will be other polarization experiments from the 
ground before the satellites fly, for example POLAR at Wisconsin University [|TT| and 



POLATRON from CalTech |12[. We illustrate the impact of the reionization signal in 
the polarization by looking at the Wisconsin experiment. It will try to measure both Q 
and U parameters with an expected sensitivity of l/xK per pixel. The instrument will 
allow measurements with a 7° FWHM at an early stage and a 1° FWHM afterwards. 
This corresponds to a gaussian window function, Wi = exp[— (/ + 0.5)^(t|], ae = 



9/2^ (2\n2) where 6 is the FWHM of the detector in radians. The predicted values 
for the Stokes parameters were calculated numerically and the spectra normalized to 
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Figure 4-4: /(/ + l)Ci/27r for both temperature (a) and polarization (b) for standard 
CDM and a model where the optical depth to recombination is Uri = 0.5 and a 
spectral index n = 1.2. 



Figure 4-5: Polarization rms fluctuations (fiK) as a function of the optical depth, Kri 
for a 7° and 1° FWHM experiments. 

COBE. The expected rms values of Q for standard CDM with no reionization are 
P(7°) = 4.8 X 10~^/iK and -P(l°) = 0.77/iK. These values, especially the large angular 
scale one, are extremely small and thus very difficult to detect. This is the reason 
why the reionized scenarios are the most promising to detect polarization. 

Reionization will not only change the polarization power spectrum but also the 
temperature one, and could suppress the acoustic peaks completely. There is some 
degree of degeneracy between the different parameters determining the CMB spectra, 
for example a reionization with a moderate optical depth will decrease the amplitude 
of the acoustic peaks but this effect may be compensated by changing the spectral 



index [|13l. In fact only an optical depth in the 0.10 — 0.20 range seems detectable from 
temperature maps alone |l|, |^, |^ . Figure [4-4| shows both polarization and temperature 
power spectra for standard CDM with a spectral index n = 1 and a reionized model 
with Kri = 0.5 but a spectral index n = 1.2. The differences in the anisotropy 
power spectra are not large, while the polarization spectra are very different. The 
rms P values in this reionized case are P{7°) = 1.2/iK and P{1°) = 1.8/iK. For the 
large angular scale experiment the difference with standard CDM is more than two 
orders of magnitude and in the one degree case is more than a factor of two. Thus 
a polarization measurement would much more easily distinguish between these two 
scenarios. 



Figure |4-5| shows the rms value of P as a function of k^, the major parameter 
determining the amplitude of the polarization perturbation. P{7°) only exceeds l/xK 
for > 0.5 but saturates quickly near 1.8/iK. On the other hand P{1°) quickly 
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raises above l/iK and reaches 3.2/iK for an optical depth of two. This means that 
even a negative detection at the 1/iK level for the one degree experiment is enough 
to rule out some models, those with optical depth Kri > 0.3. Parameters other than 
Kri do not make much difference in the height of the peaks. Table |4.1| explores the 
dependence of P{7°) and P{1°) with different cosmological parameters for a fixed 

= 1.0. Although the height of the peaks remains almost constant in these models, 
slight shifts in their location change the predicted P. The 7° rms linear polarization is 
more sensitive to the position of the first peak. The 1° experiment has the best chance 
of putting interesting constraints on a possible reionization, as the expected signal is 
greater because it is sensitive to all the power in the new peaks of the polarization 
power spectrum. A correlation analysis between the polarization in the forty pixels 
that the experiment will measure may further improve the above limits. 

In summary, the polarization of the microwave background is very sensitive to the 
ionization history of the universe and an early reionization can greatly enhance it. We 
have discussed in detail the physics behind the generation of polarization in reionized 
scenarios and the appearance of new peaks in the polarization power spectrum. We 
have identified the major parameters determining the location of these peaks, the 
ratio of distances between the observer and the reionization scattering surface to that 
between reionization and recombination. The height of the peaks is mainly function 
of Kri, the optical depth to recombination. 



4.3 The Detection of Non-Scalar Perturbations 

Primordial gravity waves produce fluctuations in the tensor component of the met- 
ric, which could result in a significant contribution to the CMB anisotropics on large 
angular scales. Unfortunately, the presence of scalar modes prevents one from clearly 
separating one contribution from another. If there are only a finite number of mul- 
tipoles where the tensor contribution is significant then there is a limit in amplitude 
beyond which tensors cannot be distinguished from random fluctuations. In a noise 
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Table 4.1: Degree of linear polarization in /iK SCDM (first row) and several cos- 
mological models all with k^i = 1.0. The value of the cosmological constant is 
such that all the above models are flat, Qtotai = 1-0. Hq is the Hubble constant 
in km sec~^Mpc~^. 



no 




Ho 


P(7°) 




1.0 


0.05 


50 


4.8110-2 


0.642 


0.7 


0.05 


50 


1.62 


2.25 


0.5 


0.05 


50 


1.67 


2.50 


0.3 


0.05 


50 


1.62 


2.25 


1.0 


0.03 


50 


1.40 


2.67 


1.0 


0.08 


50 


1.83 


2.79 


1.0 


0.10 


50 


1.91 


2.80 


1.0 


0.05 


60 


1.72 


2.79 


1.0 


0.05 


80 


1.84 


2.85 


1.0 


0.05 


100 


1.92 


2.88 



free experiment the tensor to scalar ratio T/ S needs to be larger than 0.15 to be mea- 
surable in temperature maps |[T^. Independent determination of the tensor spectral 
slope riT is even less accurate and a rejection of the consistency relation in inflation- 
ary models T/S = —Tut is only possible if \nT\ ^ (T/S)/7 [T^, [1^. Polarization 
produced by tensor modes has also been studied [|1^], but only in the small scale 
limit. In previous work correlations between the Stokes parameters Q and U have 
been used. These two variables are not the most suitable for the analysis as they 
depend on the orientation of coordinate system. In Chapter ^ we presented a full 
spherical analysis of polarization using Newman-Penrose spin-s spherical harmonic 
decomposition. We have shown that there is a particular combination of Stokes pa- 
rameters that vanishes in the case of scalar modes, which can thus be used as a probe 
of gravity waves. Here we discuss the sensitivity needed to detect this signal and com- 
pare it to the expected sensitivities of future CMB satellites We use the expressions 
in Chapters Q and ^ to evaluate the power spectra in various theoretical models. We 
use T/S as the parameter determining the amplitude of tensor polarization. Figure 
[4-6| shows the predictions for scalar and tensor contribution in standard CDM model 
with no reionization (a) and in reionized universe with optical depth of Kri = 0.2 



85 



(b). The latter value is typical in standard cosmological models |19]. We assumed 
T/S = 1 and rij' = {ug — 1) = —0.15. In the no-reionization case both tensor spectra 
peak around / ~ 100 and give comparable contributions, although the B channel is 
somewhat smaller. Comparing the scalar and tensor E channels one can see that 
scalar polarization dominates for T/S^ 1. Even though tensor contribution is larger 
than scalar at low I, the overall power there is too small to be measurable. Tensor 
reconstruction in the E channel suffers from similar drawbacks as in the case of tem- 
perature anisotropics: because of large scalar contribution, cosmic variance prevents 
one from isolating very small tensor contributions [|l^. The situation improves if 
the epoch of reionization occurred sufficiently early that a moderate optical depth to 
Thomson scattering is accumulated (Fig. [4-6|b). In this case there is an additional 



peak at low / ||20| and the relative contribution of tensor to scalar polarization in 
E channel around / = 10 is higher than around / = 100. Still, if T/S* ^ 1 cosmic 
variance again limits out ability to extract unambiguously the tensor contribution. It 
is in this limit that the importance of the B channel becomes crucial. This channel is 
not contaminated by scalar contributions and is only limited by noise, so in principle 
with sufficient noise sensitivity one can detect even very small tensor to scalar ratios. 
Moreover, a detection of signal in this channel would be a model independent detec- 
tion of non-scalar perturbations. In the following we will discuss sensitivity to gravity 
waves using both only B channel information alone and all available information. 

We can obtain an estimate of how well can tensor parameters be reconstructed by 
using only the B channel and assuming that the rest of cosmological parameters will 
be accurately determined from the temperature and E polarization measurements. 
While this test might not be the most powerful it is the least model dependent: any 
detection in B channel would imply a presence of non-scalar fluctuations and therefore 
give a signiflcant constraint on cosmological models. Because the B channel does not 
cross-correlate with either T or E ^ only its auto-correlation needs to be 
considered. A useful method to estimate parameter sensitivity for a given experiment 
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Figure 4-6: Multipole moments for the three polarization spectra for no-reionization 
case ^a ) and reionized case with optical depth of 0.2 (b). The underlying model is 
"standard CDM" with T/S = l. 



is to use the Fisher information matrix introduced in equation (|2.36|) [||, 22, 2^ 

{21 + l )/sky 

47ra2 ,„,-,,„2i-2 (dCu,\ ( dC 



(^ij = E 2 

C«,!^,<.«K,-f^gi)(^|, (4.0, 



X 



N 



where /sky is the sky coverage. Receiver noise can be parametrized by Ana'^ /N , where 
0" is the noise per pixel and is the number of pixels. Typical values are (0.15/iK)^ 
for MAP and (0.025/iK)^ for the most sensitive Planck bolometer channel in one year 
of observation. In our case the parameters Si can be T/ S and n-r, so that the matrix 
is only 2x2. The error on each parameter is given by {a'^Y^"^ if the other parameter is 
assumed to be unknown and {aa)'^^'^ if the other parameter is assumed to be known. 
Using this expression we may calculate the experiment sensitivity to these parameters. 
Current inflationary models and limits from large scale structure and COBE predict 



T / S to be less than unity. Figure |4-6| shows that the expected amplitude in this case 
is below 0.5/iK. We find that MAP is not sufficiently sensitive in the B channel to 
detect these low levels. On the other hand, Planck will be much more sensitive and 
can detect T/S > 0.2 if the tensor index ut is assumed to be known (for example 
through the consistency relation). For the underlying model with T/S = 1 one can 
determine it with an error A (T/S*) ~ 0.1. If the tensor index is not known then a 
combination of the two parameters, which corresponds to the total power under the 
B curve in Figure [4-6| , can still be determined with the same accuracy. 

Separate determination of the tensor amplitude and slope from the B channel is 
only possible in reionized models. In the no-reionization model the contribution to 
B is very narrow in / space and the leverage on rix independent of T/S is small, 
so that the correlation coefficient ai2/{aua22Y^'^ is almost always close to unity. A 
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modest amount of reionization improves the separation; in the reionized models the 
power spectrum for B is bimodal (Figure |4-6| ) and the overall signal is higher, which 
gives a better leverage on ht independent of T/S. For Uri = 0.2 the Planck errors are 
A{T/S) ~ 0.15 and Aut ~ 0.1 for the underlying model with T/S = 1. These results 
depend on the overall amplitude relative to the noise level. As long as both peaks 
can be separated from the noise one can determine the tensor slope, which allows to 
test the inflationary consistency relation. 

Combining temperature, E polarization and their cross-correlation further im- 
proves these estimates. In this case other parameters that affect scalar modes such as 
baryon density, Hubble constant or cosmological constant enter as well and the results 
become more model dependent (Chapter The Fisher information matrix has to 
be generalized to include all the parameters that can be degenerate with the tensor 
parameters. The results depend on the class of models and number of parameters 
one restricts to in the analysis, as opposed to the results based on B channel above, 
which depend only on the two main parameters that characterize the gravity wave 
production. As a typical example, for T/S = 0.1 and Uri = 0.1 one can determine 
A{T/S) = 0.05 and Arix = 0.2 with Planck 0. These errors improve further if a 
model with higher T/S or is assumed. For the same underlying model without 
using polarization the expected errors are AT/S ~ 0.26 and Arix ~ 1, significantly 
worse than with polarization. Even for MAP the limits on T/S improve by a factor 
of 2 when polarization information is included. 

To summarize the above discussion, future CMB missions are likely to reach the 
sensitivities needed to measure (or reject) a significant production of primordial grav- 
ity waves in the early universe through polarization measurements, which will vastly 
improve the limits possible from temperature measurements alone and will allow a 
test of consistency relation. The more challenging question is wether the foreground 
can be subtracted to at the required level. At low frequencies radio point sources 
and synchrotron emission from our galaxy dominate the foregrounds and both are 
polarized at a 10% level. Their contribution decreases at higher frequencies and with 
several frequency measurements one can subtract these foregrounds at frequencies 
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around 100 GHz at the required microkelvin level. At even higher frequencies dust is 
the dominant foreground, but is measured to be only a few percent polarized |^ . 



While we discussed only scalar and tensor modes, vector modes, if present be- 
fore recombination, will also contribute to both polarization channels and so could 
contaminate the signature of gravity waves. At present there are no viable cosmo- 
logical models that would produce a significant contribution of vector modes without 
a comparable amount of tensor modes. In inflationary models vector modes, even 
if produced during inflation, decay away and are not significant during recombina- 
tion. In topological defect models nonlinear sources continuously create both vector 
and tensor modes and so some of the signal in B channel could be caused by vector 
modes. Even in these models however a significant fraction of signal in B will still be 
generated by tensor modes and in any case, absence of signal in B channel would rule 
out such models. Polarization thus offers a unique way to probe cosmological models 
that is within reach of the next generation of CMB experiments. 



4.4 A Test of the Causal Structure of the Universe 

There are two competing families of models to explain the origin of the structure we 
observe in our universe: defect models, where a symmetry breaking phase transition 
generates seeds that form sub-horizon scale density fluctuations, and inflationary 
models, where a period of superluminal expansion turns quantum fluctuations into 
super-horizon density perturbations. A fundamental difference between these two 
mechanisms of structure formation is that only inflation alters the causal structure 
of the very early universe and is able to create correlations on super-horizon scales, 
while defects are causal and all correlations vanish for events where both light cones 
do not overlap. 

The COBE satellite observed correlations on angles much larger than that sub- 
tended by the horizon at decoupling {6h ~ 1.1°) in the CMB temperature. This does 
not however, imply that there were correlations on super horizon scales at decoupling 
because a time dependent gravitational potential can produce temperature fluctua- 
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tions at late-times through the integrated Sachs Wolfe effect (ISW). For example, 
cosmic string or texture models can reproduce the COBE results despite being causal 
models by generating the fluctuations at low redshift. 

Measurements of temperature fluctuations at small scales have been suggested as 
a potential test of inflation: inflationary models and most non-inflationary ones pre- 
dict different locations and relative heights for the acoustic peaks . Unfortunately, 
causality alone is insufficient to distinguish inflationary and non-inflationary temper- 
ature power spectra: causal sources that mimic exactly the inflationary pattern of 
peaks can be constructed |2^. While the predicted CMB fluctuations of the cur- 



rent family of defect models differ significantly from inflationary predictions ||2^ , it is 



useful to have model independent tests of the causal structure of the early universe. 

Polarization fluctuations are produced by Thomson scattering during the decou- 
pling of matter and radiation. Thus, unlike temperature fluctuations, they are unaf- 
fected by the ISW effect. Measurements of the polarization fluctuations are certain 
to probe the surface of last scattering. Hence, the detection of correlated polarization 
fluctuations on super-horizon scales at last scattering provides a clear signature of 
the existence of super-horizon scale fluctuations, one of the distinctive predictions of 
inflation.^ 

We will work in the initially unperturbed synchronous gauge, where the metric is 
given by ds^ = a^{T)[—dT'^ + {6ij + hij)dx'^dx^]. We will consider only perturbations 
produced by scalar modes and will solve the Einstein equations in the presence of 
sources (e.g., defects) using the stiff approximation p8[. The sources are characterized 



by their covariantly conserved stress energy tensor 6^,^. Before recombination, matter 
and radiation act as a very tightly coupled fluid, so the evolution of fluctuations can 
be described by 



Sc + -Sc = 47rG(^(l + 34)p^5,v + eoo + e) 



« N 



■^In this section we will consider the correlation function in real space (ie. as a function of the 
separation angle) rather than the usual power spectrum. By doing so we can easily express the 
causality constraint, while it would become a set of integral constraints that the power spectra have 
to satisfy in the now more usual treatment in term of Cis. 
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4 • 4 
= - -V ■ VR 

3 

VR = -(1-34)-vr--4V<5r, (4.10) 

where 9 is the trace of the spatial part of Q^u, 0/3 is the source pressure, Sr and 
vr, describe the energy density and velocity of the photon-baryon fluid and 6c is the 
energy density of cold dark matter. In synchronous gauge, the cold dark matter has 
zero velocity. The sum over is carried out over all species and cs is the sound 
speed. Temperature and polarization anisotropies seen on the sky today depend on 
6r and vr at decoupling. 

Equations ( |4.10| ) imply that the photon-baryon fluid propagates information at 
the speed of sound and thus cannot generate correlations on scales larger than the 
sound horizon. Causality on the other hand implies that the unequal time correlators 
of the sources (6^j,(r, r)6^,^(0, r')) vanishes if r > r + r'. In the absence of initial 
correlations, these two conditions together imply that (X|^^(ni)X|T-, (£12)) = if > 
29h ~ 2°, where X = 6r, vr, 9jVR and is the conformal time of decoupling. 

In the thin scattering surface approximation, equations ( [4.1U| ) are solved up to re- 



combination and then the photons free stream to the observer. The final temperature 
anisotropy in direction n on the sky is (equation |2.43| ), 



T(n) = ^1,. - n ■ vrI.. - - / dTh,,n'h^. (4.11) 



The first two terms are evaluated at the last scattering surface and the third term 
is an integral along the line of sight, the ISW effect. In non- inflationary models, the 
first two terms cannot correlate temperature fluctuations at separations larger than 
29h ~ 2° but because anisotropies can be created later through the ISW effect these 
models can have temperature correlations on larger angular scales. 
For the polarization we had (equation p^.44D , 

(Q + zf/)(n) ^ 0.17Anm'ui^diVj\r,. (4.12) 
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This equation shows that the observed polarization only depends on the state of the 
fluid at the last scattering surface. No correlations can be present in the polarization 
for separations larger than ~ 2° in non- inflationary models. An early reionization of 
the universe alters this argument, we will discuss that later. 

The correlation functions of Q and U (in their natural coordinate system) are 
given by (eq. |2.18|) , 



C^<^^^\9) = T.^-^[Cr' Fl{9) - Cl^'"^^ F^e)] (4.13) 



where ±2F;2(^,0) = ^(2/ + l)/47r [F/(0) ± F^{9)] exp (2i0). Both correlation func- 
tions receive contributions from the E and B channels. The E channel contains all 
the cosmological signal if there are no tensor or vector modes. 

We computed both C^'^'^\6) for the model proposed by Turok |2^ which has a 
clever choice of source stress energy tensor that is able to reproduce the pattern of 



peaks of inflationary standard CDM (sCDM). The results are shown in Figure [4-7 . 
We see that the inflationary model is able to produce polarization correlations on 
angular scales larger than ~ 2°, while the other model cannot. On smaller angular 



scale than shown in Figure \i-7[ the two correlation functions coincide. The difference 
between the two models is a result of the causal constraints and is insensitive to source 
evolution. It is also worth pointing out that in inflationary models the large scale 
polarization is suppressed relative to the small scale signal, so we are after a small 
effect. 

Next, we estimate the expected uncertainties in measuring bmce re- 

ceiver noise is the likely to be the dominant source of variance, we can make a simple 
estimate of the total noise: it is proportional to the number of independent pairs of 
pixels, Np, at a given separation, 6. For an experiment with a beam FWHM Ofwhm, 
Np = 1/2 X {4:71 / Oj^f^^) X {2tt9/ Ofwhm)- If o"(^q u) noise in the polarization mea- 

surement per resolution element, then the noise in the cross correlation is given by 



ACW'^) = ^2/Np afQ^u^ ^ 20 Wp^^Q^^/Ofwh^po ^here Wp^ = ctJq^^) fip„/47r. 
We can make a more accurate determination of the noise using the covariance 
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matrix of the different power spectra (equation 2.34 ): 



CoHCfj,,B)l) = ^{CiE,B)l + WpV'^'t)^ (4.14) 

which give the following variances for the correlation functions, 

+cov(cf^,^),) me)n. (4.15) 

Figure [4-7| shows the noise in each correlation, in the limit where the variances are 
dominated by receiver noise (AC^)^ = {AC'^Y and agree perfectly with our previous 
estimates. If either the cosmic variance is important or the power spectra of E and B 
differ, then the approximate estimate of the previous paragraph is not accurate and 
the full calculation should be used to estimate the noise. 

The noise in the correlation functions can be reduced by focusing on the ii^-like 
piece of the polarization. The noise in the both C^^'^\6) receives contributions 
from the variances in both E and B spectra (equation [4.15|) , but by computing both 



contributions separately we can show that the variance in E{B) makes the dominant 
contribution to AC"5(AC^). If we fiher the maps to pull out only the E component, 
then we remove not only the B signal but also some of the noise and AC^ is reduced 
almost by a factor of ~ 4. The assumption that most of the signal is in the E channel 
can be checked within the data as both E and B contributions can be measured 
separately from the maps. 

For the MAP satellite, without filtering the noise, AC^'^) 0.36/v^ ^K'^, so 
it will not be sensitive enough to detect this signal, even if we combine all of the 
three highest frequencies. However, if we filter the map to extract the E channel 
signal, then the noise in the MAP experiment drops to AC^ ~ 0.1/v^ ^^K'^■, and 
the signal should be detectable. The PLANCK satellite, with its very sensitive 
bolometers, should be able to achieve AC^'^''^^ ~ 0.003/-\/^ ^K^ and should easily be 
able to detect both and . As cosmic variance is not the dominant contribution 
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to the noise, an experiment observing a small patch of the sky could also potentially 
detect this signal. 

The temperature-polarization cross-correlation is another potential test of the 



origin of fluctuations: although ISW effects produce temperature fluctuations after 
decoupling, we still do not expect correlations between temperature and polarization 
on large angular scales for defect models. The correlations between temperature and 
polarization fluctuations directions fii and n2 are, 

{Q{h,)T{h,)) = {Q\T;) - \ rfrnlni(/^,,(r)gt), (4.16) 



T* stands for the first two terms in equation ( [4.11j ) and Q* is given by equation (|4.12| ). 



In the polarization temperature cross correlation, only the term involving the line 
of sight integration could produce correlations on large angular scales. This would 
require correlations between the late time variations of the metric and the velocity at 
last scattering. For this to occur in defect models, they must be moving very fast and 
remain coherent as they evolve from recombination to very late times. As Figure 1(c) 
shows, even Turok's causal seed model, which mimics inflation remarkable well in the 
temperature correlation does not predict any temperature-polarization correlation. 

If gravity waves, rather than scalar modes, were the dominant source of the 
anisotropics, then they could, in principle, create a cross correlations on large an- 
gular scales. However, if gravity waves were significant enough to create a large 
signal, then they would be directly detectable in the B channel. 

Figure |4-7| (c) shows the calculated values of the cross correlation together with 
the expected noise. The signal is well above the noise for MAP and Planck. The 
detection of a large angular scale cross correlation with no appreciable signal in the 
polarization B channel would put very stringent limits on the physics of models trying 
to mimic inflation. 

There is one caveat to our argument: reionization. If the universe reionized very 
early, a significant fraction of the observed polarization will come from the rescattering 
of photons at late times. In most scenarios, the fraction of rescattered photons is 
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Figure 4-7: Correlation functions for Q (a) and U (b) Stokes parameters for sCDM 
and the causal seed model discussed in the text. The noise in their determination 
is shown for both MAP and Planck. Panel (b) shows the expected noise for MAP 
if the CMB maps are filtered to include only the E channel signal. Panel (c) shows 
the cross correlation between temperature and polarization and the noise for MAP, 
the expected variance for Planck is even smaller. Each resolution element in the 
correlation function should be considered independent. 

thought to be less than ~ 20% [T^. Reionization has two effects on our argument. 



First, it reduces the amplitude of the correlation function by a factor exp(— 2^^), 
where Uri is the optical depth to decoupling (kh < 0.2). Second, it creates further 
structure in the correlation function on large angular scales. Fortunately, the effect of 
reionization can be separated from that of the primordial anisotropies: it leaves a very 
specific signature in the power spectrum, a peak at very low I that is easily distinguish 



from the I dependence expected from causality constraints alone |2y, |3y, ^ . Because 
of the form of and Ff{0), this peak produces an almost constant positive offset 

in C"^ and for angles 9 ~ 2°. Because the offset in is positive [FfiO) < 
for 9 2° and / < 70) reionization at a relatively recent epoch can never create the 
negative peak at ~ 2° predicted by inflationary models. 

There is a precise signature in C^^'^\9) on ~ 2° scales that would allow an 
unambiguous test of inflation. The signal is small, but within reach for the new gen- 
eration of experiments. The cross correlation between temperature and polarization 
is also expected to provide strong constraints that could distinguish inflation from 
non-inflationary models: this signal is much larger and will be well above the noise 
for MAP. The next generation of satellites or even polarization measurements from 
the ground could provide a definitive test of the inflationary paradigm in the relatively 
near future. 
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Chapter 5 

Predictions for Future 
Experiments 

It has long been recognized that the microwave sky is sensitive to many cosmological 
parameters, so that a high resolution map may lead to their accurate determination 
|l|, H, 1^, ^. The properties of the microwave background fluctuations are sensitive 
to the geometry of the universe, the baryon-to-photon ratio, the matter density, the 
Hubble constant, the cosmological constant, and the optical depth due to reionization 
in the universe. A stochastic background of gravitational waves also leaves an imprint 
on the CMB and their amplitude and slope may be extracted from the observations. 
In addition, massive neutrinos and a change in the slope of the primordial spectrum 
also lead to potentially observable features. 

Previous calculations trying to determine how well the various parameters could 
be constrained were based on approximate methods for computing the CMB spectra 
0. These approximations have an accuracy of several percent, which suffices for the 
analysis of the present-day data. However, the precision of the future missions will be 
so high that the use of such approximations will not be sufficient for an accurate de- 
termination of the parameters. Although high accuracy calculations are not needed at 
present to analyze the observations, they are needed to determine how accurately cos- 

^Based on M. Zaldarriaga, D. N. Spergel & U. Seljak, Astrophys. J. 488, 1 (1997) 
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mological parameters can be extracted from a given experiment. This is important 
not only for illustrative purposes but may also help to guide the experimentalists 
in the design of the detectors. One may for example address the question of how 
much improvement one can expect by increasing the angular resolution of an exper- 
iment (and by doing so increasing the risk of systematic errors) to decide whether 
this improvement is worth the additional risk. Another question of current interest 
is whether it is worth sacrificing some sensitivity in the temperature maps to gain 
additional information from the polarization of the microwave background. When 
addressing these questions, the shortcomings of approximations become particularly 
problematic. The sensitivity to a certain parameter depends on the shape of the 
likelihood function around the maximum, which in the simplest approach used so far 
is calculated by differentiating the spectrum with respect to the relevant parameter. 
This differentiation strongly amplifies any numerical inaccuracies and almost always 
leads to an unphysical breaking of degeneracies among parameters and misleadingly 
optimistic results. 

Previous analysis of CMB sensitivity to cosmological parameters used only tem- 
perature information. However, CMB experiments can measure not just the tempera- 
ture fluctuations, but also even weaker variations in the polarization of the microwave 
sky. Instead of one power spectrum, one can measure up to four and so increase the 
amount of information in the two-point correlators ^ |^. Polarization can provide 
particularly useful information regarding the ionization history of the universe and 
the presence of a tensor contribution p, |^. Because these parameters are partially 
degenerate with others, any improvement in their determination leads to a better 
reconstruction of other parameters as well. The two proposed satellite missions are 
currently investigating the possibility of adding or improving their ability to measure 
polarization, so it is particularly interesting to address the question of improvement 
in the parameter estimation that results from polarization. 

The purpose of this Chapter is to re-examine the determination of cosmological 
parameters by CMB experiments in light of the issues raised above. It is particularly 
timely to perform such an analysis now, when the satellite mission parameters are 
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roughly defined. We use the best current mission parameters in the calculations 
and hope that our study provides a useful guide for mission optimization. As in 
previous work , we use the Fisher information matrix to answer the question of how 
accurately parameters can be extracted from the CMB data. This approach requires 
a fast and accurate method for calculating the spectra and we use the CMBFAST 
code developed in Chapter ^ |jlO| with an accuracy of about 1%. We test the Fisher 
information method by performing a more general exploration of the shape of the 
likelihood function around its maximum and find that this method is sufficiently 
accurate for the present purpose. 



The outline of this chapter is the following: in Section we present the methods 
used, reviewing the calculation of theoretical spectra and the statistical methods 



to address the question of sensitivity to cosmological parameters. In Section 5.2 



we investigate the parameter sensitivity that could be obtained using temperature 
information only and in Section we repeat this analysis using both temperature 
and polarization information. In Section ^.4| , we explore the accuracy of the Fisher 
method by performing a more general type of analysis and investigate the effects of 
prior information in the accuracy of the reconstruction. We present our conclusions 



in Section 5.5 



5.1 Methods 

In this section, we review the methods used to calculate the constraints on different 
cosmological parameters that could be obtained by the future CMB satellite experi- 
ments. We discuss the Fisher information matrix approach, as well as the more general 
method of exploring the shape of the likelihood function around the minimum. 



5.1.1 The Fisher information matrix 

The Fisher information matrix is a measure of the width and shape of the likelihood 
function around its maximum. Its elements are defined as expectation values of 
the second derivative of a logarithm of the likelihood function with respect to the 
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corresponding pair of parameters. It can be used to estimate the accuracy with 
which the parameters in the cosmological model could be reconstructed using the 
CMB data 0, [ll| . If only temperature information is given then for each I a derivative 



of the temperature spectrum Cti with respect to the parameter under consideration 
is computed and this information is then summed over all I weighted by Cov^^ (C*|.J. 
In the more general case implemented here, we have a vector of four derivatives and 
the weighting is given by the inverse of the covariance matrix, 



a. 



s:s:l^cov-(&„cv,)^. (5.1) 

I X,Y 



Here aij is the Fisher information matrix, Cov~^ is the inverse of the covariance 
matrix, Si are the cosmological parameters one would like to estimate and X, Y stands 
for T,E,B,C. For each /, one has to invert the covariance matrix and sum over X 
and Y. The derivatives were calculated by finite differences and the step was usually 
taken to be about 5% of the value of each parameter. We explored the dependence 
of our results on this choice and found that the dependence is less than 10%. This 
indicates that the likelihood surface is approximately Gaussian, an assumption of the 
Fisher matrix method that only consideres the curvature around the maximum of the 
likelihood. Further tests of this assumption are discussed in Section B^. 



The full covariance matrix between the power spectrum estimators was presented 
in p, H, |[. The diagonal terms are given by 



[ZL + I) J sky 

Cov(C'|,) = {CEi + Wp'Bry 

[ZL + I) J sky 

CoviCl,) = ] {CBi + Wp'Bf'f 

[ZL + l)Jsky 

Cov(C'5,) = ^—^[cI, + {Cti + w^^By'){Cei + Wp'BY%, (5.2 

[ZL + I) J sky '- 
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while the non-zero off diagonal terms are 



(21 + 1)/, 



ky 



2 

'CI 



Com{CtiCci) = Cci{CTi + WT^Br'') 

[ZL + I) J sky 

Cov{CeiCci) = i^.^,,r CciiCEi + Wp'Br'). (5.3) 

[ZL + l)Jsky 

We have defined W^r^p^ = Ana'^j, p-^/Npi^ where and ap are noise per pixel in the 
temperature and either Q or U polarization measurements (they are assumed equal) 
and Npij. is the number of pixels. We will also assume that noise is uncorrelated be- 
tween different pixels and between different polarization components Q and U. This 
is only the simplest possible choice and more complicated noise correlations arise if all 
the components are obtained from a single set of observations. If both temperature 
and polarization are obtained from the same experiment by adding and differentiating 
the two polarization states, then = ap/2 and noise in the temperature is uncorre- 
lated with the noise in polarization components. The window function Bf^ accounts 
for the beam smearing and in the Gaussian approximation is given by Bf^ = exp /^cr^, 
with cTfo measuring the width of the beam. We introduced fsky as the fraction of the 
sky that can be used in the analysis. In this chapter we assume fsky = 0.8. It should 
be noted that equations ( |5.2| ) and ( ^.3| ) are valid only in the limit of uniform sky 
coverage. 

Both satellite missions will measure in several frequency channels with different 
angular resolutions: we combine them using W(^t,p) = J2'>^Itp)^ where subscript c 
refers to each channel component. For the MAP mission we adopt a noise level w^^ = 
(O.ll/ii^)^ and Wp^ = (O.lS/ii^)^ for the combined noise of the three highest frequency 
channels, with conservatively updated MAP beam sizes: 0.53°, 0.35° and 0.25°. These 
beam sizes are smaller than those in the MAP proposal and represent improved 
estimates of MAP's resolution. The most recent estimates of MAP's beam sizes are 
even smalleif): 0.47°, 0.35° and 0.21°. For Planck, we assume w^^ = {O.OllfiKy 



^See the MAP homepage at http:/ /map. gsfc.nasa.gov. 
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and Wp^ = (0.025/xii')^, and combining 140 GHz and 210 GHz bolometer channels. 
For Planck's polarization sensitivity, we assume a proposed design in which eight 
out of twelve receivers in each channel have polarizers. The angular resolution at 
these frequencies is 0.16° and 0.12° FWHM. We also explore the possible science 
return from an enhanced bolometer system that achieves polarization sensitivity of 
Wp^ = {O.Q15fiKy. 

In our analysis, we are assuming that foregrounds can be subtracted from the 
data to the required accuracy. Previous studies of temperature anisotropies have 
shown that this is not an overly optimistic assumption at least on large angular scales 



(e.g. ||T2|]). On smaller scales, point source removal as well as secondary processes 
may make extracting the signal more problematic. This would mostly affect our 
results on Planck, which has enough angular resolution to measure features in the 
spectrum to / ~ 3000. For this reason, we compared the results by changing the 
maximum / from 3000 to 1500. We find that they change by less than 30%, so that 
the conclusions we find should be quite robust. Foregrounds for polarization have 
not been studied in detail yet. Given that there are fewer foreground sources of 
polarization and that polarization fractions in CMB and foregrounds are comparable, 
we will make the optimistic assumption that the foregrounds can be subtracted from 
the polarization data with sufficient accuracy as well. However, as we will show, 
most of the additional information from polarization comes from very large angular 
scales, where the predicted signal is very small. Thus, one should take our results 
on polarization as preliminary, until a careful analysis of foreground subtraction in 
polarization shows at what level can polarization signal be extracted. 

The inverse of the Fisher matrix, a~^, is an estimate of the covariance matrix 



between parameters and y approximates the standard error in the estimate of 

the parameter Sj. This is the lower limit because Cramer- Rao inequality guarantees 
that for an unbiased estimator the variance on the i-th parameter has to be equal to 



or larger than \J{a~^)ii. In addition to the diagonal elements of a"^, we will also use 
2x2 submatrices of to analyze the covariance between various pairs of parameters. 
The Fisher matrix depends not only on the experiment under consideration, hut also 
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on the assumed family of models and on the number of parameters that are being 
extracted from the data. To highlight this dependence and to assess how the errors 
on the parameters depend on these choices we will vary their number and consider 
several different underlying models. 

5.1.2 Minimization 

The Fisher information matrix approach assumes that the shape of the likelihood 
function around the maximum can be approximated by a Gaussian. In this section, we 
drop this assumption and explore directly the shape of the likelihood function. We use 
the PORT optimization routines to explore one direction in parameter space at a 
time by fixing one parameter to a given value and allowing the minimization routine 
to explore the rest of parameter space to find the minimum of = Y^i J2x,y(^^i "~ 
Cxi)Cov^^{CxiCyi){Cyi — Cyi), where Cxi denotes the underlying spectrum. The 
value of a function of this parameter can be compared directly with the Fisher 
matrix prediction, Ax^ = {si — s*Y / {a^^)ii, where s* is the value of the parameter 
in the underlying model. This comparison tests not only the shape of the likelihood 
function around the maximum but also the numerical inaccuracies resulting from 
differentiating the spectrum with respect to the relevant parameter. The minimization 
method is also useful for finding explicit examples of degenerate models, models with 
different underlying parameters but almost indistinguishable spectra. 

The additional advantage of the minimization approach is that one can easily 
impose various prior information on the data in the form of constraints or inequali- 
ties. Some of these priors may reflect theoretical prejudice on the part of the person 
performing the analysis, while others are likely to be less controversial, such as the 
requirement that matter density, baryon density and optical depth are all positive. 
One might also be interested in incorporating priors into the estimation to take other 
astrophysical information into account, e.g. the limits on the Hubble constant or Q^n 
from the local measurements. Such additional information can help to break some 
of the degeneracies present in the CMB data, as discussed in Section 3. Note that 
prior information on the parameters can also be incorporated into the Fisher matrix 
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analysis, but in its simplest formulation only in the form of gaussian constraints and 
not in the form of inequalities. 

The main disadvantage of this more general analysis is the computational cost. 
At each step the minimization routine has to compute derivatives with respect to all 
the parameters to find the direction in parameter space towards the minimum. If the 
initial model is sufficiently close to the minimum, then the code typically requires 
5-10 steps to find it and to sample the hkelihood shape this has to be repeated for 
several values of the parameter in question (and also for several parameters). This 
computational cost is significantly higher than in the Fisher matrix approach, where 
the derivatives with respect to each parameter need to be computed only once. 

5.2 Constraints from temperature data 

In this section, we investigate how measurements of the CMB temperature anisotropics 
alone can constrain different cosmological parameters. The models studied here are 
approximately normahzed to COBE, which sets the level of signal to noise for a 
given experiment. We will start with models in a six-dimensional parameter space 
Sq = {C2^\h,^lx,^b, '^rijiT's), where the parameters are respectively, the amplitude 
of the power spectrum for scalar perturbations at / = 2 in units of /jK^, the Hubble 
constant in units of 100 km/s/Mpc, the cosmological constant and baryon density 
in units of critical density, the reionization optical depth and the slope of primor- 
dial density spectrum. In models with a non-zero optical depth, we assume that the 
universe is instantaneously and fully reionized, so that the ionization fraction is 
before the redshift of reionization Zri and 1 afterwards. We limit to this simple case 
because only the total optical depth be usefully constrained without polar- 

ization information. We discuss the more general case when we discuss polarization 
below. 

The underlying model is standard CDM with Sq = (796, 0.5, 0, 0.05, 0.05, 1.0). Our 
base model has an optical depth of 0.05, corresponding to the epoch of reionization 
at Zri ~ 13. Models include gravity waves, fixing the tensor amplitude using the con- 
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sistency relation predicted by inflation T/S = — Trir and assuming a relation between 
the scalar and tensor spectral slopes = — 1 for < 1 and = otherwise, 
which is predicted by the simplest models of inflation. The results for MAP are sum- 
marized in Table It is important to keep in mind that the parameters are highly 
correlated. By investigating confidence contour plots in planes across the parameter 
space, one can identify combinations of parameters that can be more accurately de- 
termined. Previous analytical work [jl4|, [15| showed that the physics of the acoustic 
oscillations is mainly determined by two parameters, Q^h'^ and VLmh'^-, where VLm is the 
density of matter in units of critical density. There is an approximately flat direction 
in the three-dimensional space of Vtb and Vim'- for example, one can change VLm and 
adjust h and Vtb to keep and Qmh"^ constant, which will not change the pattern of 
acoustic oscillations. This degeneracy can be broken in two ways. On large scales, the 
decay of the gravitational potential at late times in Qm 7^ 1 models (the so called late 
time integrated Sachs- Wolfe or ISW term) produces an additional component in the 
microwave anisotropy power spectrum, which depends only on Q\ Because the 
cosmic variance (finite number of independent multipole moments) is large for small 
/, this effect cannot completely break the degeneracy. The second way is through the 
change in the angular size of the acoustic horizon at recombination, which shifts all 
the features in the spectrum by a multiplicative factor. Around = 1? this shift is 
a rather weak function of Qm and scales approximately as leading to almost no 

effect at low /, but is increasingly more important towards higher /. MAP is sensitive 
to multipole moments up to / ~ 800, where this effect is small. Consequently MAP's 
ability to determine the cosmological constant will mostly come from large scales and 
thus will be limited by the large cosmic variance. Planck has a higher angular res- 
olution and significantly lower noise, so it is sensitive to the change in the angular 
size of the horizon. Because of this Planck can break the parameter degeneracy and 
determine the cosmological constant to a high precision, as shown in Table ^.2| . 

Figure |5-l| a shows the confidence contours in the Qm — h and — h planes. 
The error ellipses are significantly elongated along the lines AQh/Qh + 2.1 Ah /h = 
and AQm/^m + 3.0Ah/h = 0. The combinations Q^h"^'^ and ^Imh^ are thus better 
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Figure 5-1: MAP confidence contours (68% and 95%) for models in tlie six param- 
eter space (a) and seven parameter space witli T/S added as a free parameter (b). 
Parameters are normalized to their value in the underlying model denoted with an 
asterisk. 

determined than the parameters fi^, fi;, and h themselves, both to about 3% for MAP. 
It is interesting to note that it is Vtmh^ rather than Qmh'^ that is most accurately 
determined, which reflects the fact that ISW tends to break the degeneracy discussed 
above. However, because the ISW effect itself can be mimicked by a tilt in the 
spectral index the degeneracy remains, but is shifted to a different combination of 
parameters. One sigma standard errors on the two physically motivated parameters 
are /S.{yt}jh^) / Vt^^h^ 3% and A{Qrnh^)/^mh^ ~ 5%. The fact that there is a certain 
degree of degeneracy between the parameters has already been noted in previous work 
(e.g. i). 

Another approximate degeneracy present in the temperature spectra is between 

(S) 

the reionization optical depth and amplitude C2 ■ Reionization uniformly sup- 
presses the anisotropics from recombination by e~'^'''. On large angular scales, new 
anisotropics are generated during reionization by the modes that have not yet entered 
the horizon. The new anisotropics compensate the e"'''"' suppression, so that there 
is no suppression of anisotropics on COBE scales. On small scales, the modes that 
have entered the horizon have wavelengths small compared to the width of the new 
visibility function and so are suppressed because of cancellations between positive 
and negative contributions along the line of sight and become negligible. The net 

result is that on small scales the spectrum is suppressed by e~'^'^'"^ compared to the 

(s) 

large scales. To break the degeneracy between C2 and Kri one has to be able to 
measure the amplitude of the anisotropics on both large and small scales and this is 
again limited on large scales by cosmic variance. Hence one cannot accurately deter- 
mine the two parameters separately, while their combination C2e~^'*'' is much better 



constrained. Figure 5-2 shows that indeed the error ellipsoid is very elongated in the 



direction AC2/C2 — O-lAnri/nri = 0, which corresponds to the above combination for 
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Figure 5-2: Confidence contours (68% & 95%) in the — plane for models in 
the six parameter space described in the text with (dotted lines) or without (solid 
lines) polarization information. 

Uri = 0.05. 

We now allow for one more free parameter, the ratio of the tensor to scalar 
quadrupole anisotropy T/S, fixing the tensor spectral index tit using the consistency 
relation predicted by inflation T/S = —Iht but not assuming a relation between ht 



and Us- The variances for MAP are again summarized in Table A comparison 
with the previous case shows that most variances have increased. The error bars for h 
and are approximately five times larger than before while that for has increased 
by a factor of six and that for Vt^h'^ by almost four. On the other hand, the error 
bar for remains unchanged. It is instructive to look again at the contour plots in 



the f2m ~ h and ^Ib — h planes shown in Figure |5-l| b. The degeneracy on individual 
parameters is significantly worse because the large angular scale amplitude can now 
be adjusted freely with the new extra degree of freedom, the tensor to scalar ratio 
T/S. This can therefore compensate any large scale ISW term and so the degeneracy 
between h, Q\ and Qi, cannot be broken as easily. However, a combination of the two 
parameters is still well constrained as shown in Figure |5- 1| . The degenerate lines are 
now given by AQ^/Qh + 1.66Ah/h = and AQm/^m + 3.0Ah/h = 0, with relative 
errors A{^lbh^-'^'^)/^lbh^-*^'^ ^ 4% and A{Vtmh^) / Vtmh^ ^ 4%, almost unchanged from 
the 6-parameter case. On the other hand for the physically relevant parameters VLbh"^ 
and VLmh"^ we now have AiVlbh'^) / Vt^h'^ 10% and AiVlmh'^) /^mh^ ~ 25%, which is 
worse than before. This example indicates how the errors on individual parameters 
can change dramatically as we add more parameters while certain combinations of 
them remain almost unaffected. 

The output of a minimization run trying to fit sCDM temperature power spectra 
with models constrained to have = 0.6 shows how different parameters can be 
adjusted in order to keep the power spectrum nearly the same. The minimization 
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Figure 5-3: Power spectra of (a) temperature and (b) polarization for two models 
that will be degenerate for MAP if only temperature information is used. The model 
with f^A = 0.6 is the result of the minimization relative to the sCDM for models 
constrained to have f^A = 0.6. Polarization helps to break this degeneracy. 

program found the model sj = (610, 0.67, 0.6, 0.03, 0.09, 1.1, 0.68) where the last num- 
ber now corresponds to the T/S ratio, as a model almost indistinguishable from the 
underlying one. The two models differ by Ax^ = 1.8 and are shown in Figure |5-3| . 

It is interesting to analyze how different parameters are adjusted to reproduce 
the underlying model. By adding gravity waves and increasing both the spectral 
index and the optical depth, the ISW effect from the cosmological constant can be 
compensated so that it is only noticeable for the first couple of C;'s. The relatively 
high amount of tensors {T/S ~ 0.7) lowers the scalar normalization and thus the 
height of the acoustic peaks, which is compensated by the increase in the spectral 
index to = 1.1 and the decrease of Qmh'^ from 0.25 to 0.18. The latter moves 
the matter radiation equality closer to recombination increasing the height of the 
peaks. This is the reason why the degeneracy line is not that of constant fim^^ as 
Figure |5-1| shows. Changes in Qmh"^ change the structure of the peaks and this can 
be compensated by changing other parameters like the optical depth or the slope of 
the primordial spectrum. This cannot be achieved across all the spectrum so one 
can expect that the degeneracy will be lifted as one increases the angular resolution, 
which is what happens if Planck specifications are used (Table p.lj ). 

Note that the amount of gravity waves introduced to find the best fit does not 
follow the relation between tit and Ug predicted by the simplest inflationary models 
discussed previously: for = 1.1 no gravity waves are predicted. This explains 
why the addition of T/S as a free parameter increases the sizes of most error bars 
compared to the 6-parameter case. 

While the two models shown in Figure have very similar temperature anisotropy 



spectra, they make very different astronomical predictions. Figure |5-4| shows the mat 



ter power spectra of the two models. An interesting effect is that the two models are 
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Figure 5-4: Hubble diagram for Type la supernovae (a) and CDM linear power spectra 
(b) for sCDM and the = 0.6 model described in the text. 



nearly identical on the scale of = O.lh Mpc~^, which corresponds to / ~ /ctq ~ 600, 
the / range where MAP is very sensitive and gravity waves are unimportant. How- 
ever, the two models differ significantly on the O.Ol/i Mpc^^ scale and the power 
spectrum shape is very different. The current observational situation is still contro- 



versial (e.g. 113)5 but measurements of the spectrum by the Sloan Digital Sky Survey 
(SDSS) should significantly improve the power spectrum determination. The models 
also make different predictions for cluster abundances: the matter dominated model 
has o"8^m^ = 1-2, while the vacuum dominated model has CgfiJ^^ = 0.8. Analysis of 
cluster X-ray temperature and luminosity functions suggests crgfi^^^ = 0.5 ± 0.1 [|18] 



inconsistent with both of the models in the figure. These kind of measurements can 
break some of the degeneracies in the CMB data. 

Observations of Type la supernovae at redshifts z ~ 0.3 — 0.6 is another very 
promising way of measuring cosmological parameters. This test complements the 
CMB constraints because the combination of Qm and Q\ that leaves the luminosity 
distance to a redshift z ~ 0.3 — 0.6 unchanged differs from the one that leaves the 
position of the acoustic peaks unchanged. Roughly speaking, the SNe observations 
are sensitive to go — — ^a, while the CMB observations are sensitive to the 

luminosity distance which depends on a roughly orthogonal combination, Q^n + ^a- 
Figure |5-4|a shows the apparent magnitude vs. redshift plot for supernovae in the 



two models of Figure |5-2| . The analysis in |[T£] of the first seven supernovae already 
excludes the Qa = 0.6 model with a high confidence. However, it remains to be 
seen however whether this test will be free of systematics such as evolutionary effects 
that have plagued other classical cosmological tests based on the luminosity-redshift 
relation. 

Finally, we may also relax the relation between tensor spectral index and its 
amplitude, thereby testing the consistency relation of inflation. For MAP, we studied 
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two sCDM models, one with T/S = 0.28 and one with T/S = 0.1, but with Kri = 0.1, 
for Planck we only used the latter model. Table |5.1| summarizes the obtained one 
sigma limits. A comparison between the T/S = 0.28 model and previous results for 
sCDM with seven parameters shows that the addition of ht as a new parameter does 
not significantly change the expected sensitivities to most parameters. The largest 
change, as expected, is for the tensor to scalar ratio. We now find AT/S ~ 0.7 which 
means that the consistency relation will only be poorly tested from the temperature 
measurements. If T/S = O.landKri = 0.1 most error bars are smaller than if T/S = 
0.28andK,ri = 0.05 case. The reason for this is that the higher value of the optical 
depth in the underlying model makes its detection easier and this translates to smaller 
error bars on the other parameters. The only exception is ny, which has a significantly 
higher error if T/S = 0.1 than if T/S = 0.28 as expected on the basis of the smaller 
contribution of tensor modes to the total anisotropics. A comparison between the 
expected MAP and Planck performances for the T/ S = 0.1 model shows that Planck 
error bars are significantly smaller. For h, Q^h'^ and the improvement is by a factor 
of 10 — 20, while for T/ S and C2^^ by a factor of 2 — 3. The limits on and n-r remain 
nearly unchanged, reflecting the fact that these parameters are mostly constrained on 
large angular scales which are cosmic variance and not noise/resolution limited. It is 
for these parameters that polarization information helps significantly, as discussed in 
the next section. 

The accuracy with which certain parameters can be determined depends not only 
on the number of parameters but also on their "true" value. We tested the sensitivity 
of the results by repeating the analysis around a cosmological constant model Sg = 
(922, 0.65, 0.7, 0.06, 0.1, 1.0, 0.1, 0.0), where the last number corresponds to the tensor 
spectral index n^. Results for MAP specifications are given in Table ^.1| . The most 
dramatic change is for the cosmological constant, which is a factor of ten better 
constrained in this case. This is because the underlying model has a large ISW effect 
which increases the anisotropics at small /. This cannot be mimicked by adjusting the 
tensors, optical depth and scalar slope as can be done if the slope of the underlying 



model is flat, such as for sCDM model in Figure |5-3| . Because of the degeneracy 
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between f^A and h , a better constraint on the former will also improve the latter, as 
shown in Table |5.1| . Similarly, because a change in Q\ affects T/ S, Uri and on large 
scales, the limits on these parameters will also change. On the other hand, errors on 
Qbh"^, C'^^ and ut do not significantly change. This example clearly shows that the 
effects of the underlying model can be rather significant for certain parameters, so one 
has to be careful in quoting the numbers without specifying the "true" parameters of 
the underlying model as well. 

So far we only discussed fiat cosmological models. CMBFAST can compute 
open cosmological models as well, and we will now address the question of how 
well can curvature be determined using temperature data. We consider models 
in a six parameter space = (C2 , h,Q,Qf,, Kri,ns), with no gravity waves and 
where Qa = 0, so that Q = Q^n- We will consider as the underlying model = 



(1122, 0.65, 0.4, 0.06, 0.05, 1.0). Fisher matrix results are displayed in Table |5Tl| . Within 
this family of models Q can be determined very precisely by both MAP and Planck 
due to its effect on the position of the acoustic peaks. This conclusion changes drasti- 
cally if we also allow cosmological constant, in which case Q = Qm + ^a- Both Q and 
Qa change the angular size of the sound horizon at recombination so it is possible to 
change the two parameters without changing the angular size, hence the two param- 
eters will be nearly degenerate in general. We will discuss this degeneracy in greater 
detail in the next section, but we can already say that including both parameters in 
the analysis increases the error bar on the curvature dramatically. 

To summarize our results so far, keeping in mind that the precise numbers depend 
on the underlying model and the number of parameters being extracted, we may rea- 
sonably expect that using temperature information only MAP (Planck) will be able 
to achieve accuracies of ACf ^ ~ 0.5(0.1), Ah ~ 0.1(0.006), A^a ~ 0.6(0.03), 
A{nbh^) /Qbh'^ ~ 0.1(0.008), Ak^ ~ 0.1(0.1), An, ~ 0.07(0.006), A(T/5) ~ 0.7(0.3) 
and AriT ~ 1(1)- It is also worth emphasizing that there are combinations of 
the parameters that are very well constrained, e.g. A{Qrnh^) /^mh^ ~ 0.04 and 
A{c!f^)/C2^^ — 2AKri ~ 0.05. For the family of models with curvature but no cos- 
mological constant, MAP (Planck) will be able to achieve AQ ~ 0.007(0.0006), de- 
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termining the curvature of the universe with an impressive accuracy. 

These results agree quahtatively, but not quantitatively, with those in The 

(S) 

discrepancy is most significant for C2 \ h and ^Ia, for which the error bars obtained 
here are significantly larger. The limit we obtain for Q^h"^ is several times smaller than 
that in [Q, while for the rest of the parameters the results agree. The use of different 
codes for computing model predictions is probably the main cause of discrepancies 
and emphasizes the need to use high accuracy computational codes when performing 
this type of analysis. 

5.3 Constraints from temperature and polariza- 
tion data 

In this section, we consider the constraints on cosmological parameters that could 
be obtained when both temperature and polarization data are used. To generate 
polarization, two conditions have to be satisfied: photons need to scatter (Thom- 
son scattering has a polarization dependent scattering cross-section) and the angular 
distribution of the photon temperature must have a non-zero quadrupole moment. 
Tight coupling between photons and electrons prior to recombination makes the pho- 
ton temperature distribution nearly isotropic and the generated polarization very 
small, specially on scales larger than the width of the last scattering surface. For this 
reason polarization has not been considered previously as being important for the 
determination of cosmological parameters. However, early reionization increases the 
polarization amplitude on large angular scales in a way which cannot be mimicked 
with variations in other parameters 0. The reason for this is that after recombi- 
nation the quadrupole moment starts to grow due to the photon free streaming. If 
there is an early reionization with sufficient optical depth, then the new scatterings 
can transform this angular anisotropy into polarization. This effect dominates on the 
angular scale of the horizon when reionization occurs. It will produce a peak in the 
polarization power spectrum with an amplitude proportional to the optical depth, 
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Figure 5-5: Confidence contours (68% & 95%) in the (a) Qb — h plane and (b) Qm — h 
plane for models in the seven parameter space described in the text with or without 
using polarization information. 



Kri, and a position / ~ 2y/zri, where Zri is the redshift of reionization. 

We first consider the six parameter space described in the previous section. Table 
|5.2| contains the one sigma errors on the parameters for MAP specifications. Com- 
pared to the temperature case, the errors improve particularly on the amplitude, the 

reionization optical depth and the spectral index n^. Figure |5-2| shows the confi- 

( S) 

dence contours for C2 and Kri with and without polarization. One can see from 
this figure how the information in the polarization breaks the degeneracy between 
the two parameters by reducing the error on Uri, but does not really improve their 
non-degenerate combination, which is well determined from the temperature data 
alone. 

We now allow for one more parameter, T/S. Again, polarization improves the 
errors on most of the parameters by a factor of two compared to the no-polarization 
case, as summarized in Table |5.2| . The optical depth and the amplitude are better 
constrained for the same reason as for the six parameter model discussed above. 
Without polarization, the extra freedom allowed by the gravity waves made it possible 
to compensate the changes on large angular scales caused by the ISW, while the 
amplitude of small scale fluctuations could be adjusted by changing the optical depth 
and the spectral index. Changing Ug also changes the slope on large angular scales, 
compensating for the change caused by the ISW. When polarization is included, a 
change on the optical depth produces a large effect in the spectrum: see the model 
with = 0.6 in Figure |5-2|, which has Kri = 0.1. The difference in between the 



two models in Figure |5-2| becomes 10 instead of 1.8 (for MAP). 

Figure |5-5| shows how the confidence contours in the ~ h and Q^ — h planes are 
improved by including polarization. The 95% confidence contour corresponds roughly 
to the 68% confidence contour that could be obtained from temperature information 
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alone, while the orientation of the error elhpsoids does not change. As before the 
well determined combination is constrained from the temperature data alone. The 
constraints on tensor parameters also improve when polarization is included. Again, 
this results from the better sensitivity to the ionization history, which is partially 
degenerate with the tensor contribution, as discussed in the previous section. The B 
channel to which only gravity waves contribute is not providing additional information 
in the model with T/S = 0.28 for MAP noise levels. Even in a model with T/S = 1 
the B channel does not provide additional information in the case of MAP. 

With its very sensitive bolometers, Planck has the potential to detect the B chan- 
nel polarization produced by tensor modes: the B channel provides a signature free of 
scalar mode contribution ||^, |^ . However, it is important to realize that even though 
for a model with T/S = 1 only 20% of the sensitivity of Planck to tensor modes is 
coming from the B channel. Planck can detected primordial gravity waves in models 
with T/S ^ 0.3 in the B channel alone. If the bolometer sensitivities are improved 
so that Wp^ = (0.015/ifr)^, then Planck can detect gravity waves in the B channel 
even if T/S ~ 0.1. We also analyzed the 8-parameter models presented in Table 
For the models with T/S = 0.28 and T/S = 0.1, MAP will not have sufficient 
sensitivity to test the inflationary consistency relation T/S = —Ttit- Planck should 
have sufficient sensitivity to determine rijn with an error of 0.2 if T/S ~ 0.1, which 
would allow a reasonable test of the consistency relation. 

Polarization is helping to constrain most of the parameters mainly by better 
constraining K^i and thus removing some of its degeneracies with other parameters. 
Planck will be able to determine not only the total optical depth through the ampli- 
tude of the reionization peak but also the ionization fraction, Xe, through its position. 
To investigate this, we assumed that the universe reionized instantaneously at Zri 
and that Xe remains constant but different from 1 for z < Zj..-/). The results given in 
Table |5.2| indicate that be determined with an accuracy of 15%. This together 

with the optical depth will be an important test of galaxy formation models which 
at the moment are consistent with wildly different ionization histories and cannot be 



probed otherwise ||2^, We also investigated the modified Planck design, where 
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Figure 5-6: Relative improvement in the parameter estimation as a function of the 
maximum / up to which polarization information is used for the MAP mission. 

Figure 5-7: Same as figure 6 but for Planck mission parameters. 

both polarization states in bolometers are measured. An improved polarization noise 
of Wp^ = {O.OlbfiKy for Planck will shrink the error bars presented in Table ^.2| by 
an additional 6 — 20 %. Error bars on Q\ and Qbh"^ are reduced by 20%, those in h, 
Kri and Xe 10 — 15 % and for T/S, tit and Us the improvement is approximately 6%. 

We can examine in more detail how polarization helps to constrain different cosmo- 
logical parameters by investigating the angular scales in the polarization power spec- 
tra that contribute the most information. To do so, we will consider the T/S* = 0.1 
model and perform a Fisher matrix analysis that includes all the temperature infor- 
mation, but polarization information only up to maximum /. Figures |5-6| (MAP) and 
[5-7| (Planck) show the increase in accuracy as a function of maximum / for various 
parameters. In the case of Planck, we added the ionization fraction after reionization 
as another parameter. Most of the increase in information is coming from the low / 
portion of polarization spectrum, primarily from the peak produced by reionization 
around / ~ 10. The first acoustic peak in the polarization spectra at / ~ 100 explains 
the second increase in information in the MAP case. The better noise properties and 
resolution of Planck help to reach the higher / polarization acoustic peaks, which add 
additional information for constraining h, Qi,h'^ and ^a- For Planck, on the other 
hand, some of the degeneracies will already be lifted in the temperature data alone 
and so less is gained when polarization data is used to constrain the ionization history. 

An interesting question that we can address with the methods developed here 
is to what extent is one willing to sacrifice the sensitivity in temperature to gain 
sensitivity in polarization. A specific example is the 140 GHz channel in Planck, 
where the current proposal is to have four bolometers with no polarization sensitivity 
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and eight bolometers which are polarization sensitive so that they transmit only one 
polarization state while the other is being thrown away. One can compare the results 
of the Fisher information matrix analysis for this case with the one where all twelve 
detectors have only temperature sensitivity, but with better overall noise because no 
photons are being thrown away. The results in this case for the 8 parameter model 
with T/S = 0.1 are 10-20% better than the results given in the fifth column of Table 
]l|. These results should be compared to the same case with polarization in Table 
The latter case is clearly better for all the parameters, especially for those that 
are degenerate with reionization parameters, where the improvement can be quite 
dramatic. Based on this example it seems clear that it is worth including polarization 
sensitivity in the bolometer detectors, even at the expense of some sensitivity in the 
temperature. However, it remains to be seen whether such small levels of polarization 
can be separated from the foregrounds. 

The Fisher matrix results for the six parameter open models are presented in 



Table |5.2| . As expected polarization improved the constraints on C2 and the 
most. So far we have explicitly left out of the analysis; as discussed in Section 3 
the positions of the peaks depends on both Q and Q/y and it is possible to change the 
two parameters without changing the spectrum. For any given value of Qm we may 
adjust h and fi^ to keep Q^h"^ and VLmh?' constant, so that acoustic oscillations will not 
change. If we then in addition adjust also VL\ to match the angular size of the acoustic 
features, then the power spectra for two models with different underlying parameters 
remain almost unchanged. As mentioned in previous section the effect of on the 
positions of the peaks is rather weak around Vt^n = 1 and the peak positions are mostly 
sensitive to the curvature Q. The lines of constant Ipeak, the inverse of the angular 
size of acoustic horizon, roughly coincide with those of constant Q near flat models, 
making it possible to weigh the universe using the position of the peaks. In the more 
general case, it is not Q that can be determined from the CMB observations, but a 
particular combination of Qm and Q\ that leaves Ipeak unchanged. Figure |^ shows 
confidence contours in the Qrn ~ plane. The contours approximately agree with 
the constant Ipeak (dotted) line, which around Vt^ = 1 coincides with the constant Q 
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Figure 5-8: Confidence contours (68% & 95%) in the — plane for open models 
in the seven parameter space described in the text. The dots show the positions of 
the minima found by the minimization routine when constrained to a subspace of 
constant VL. 



Figure 5-9: Power spectra for (a) temperature and (b) polarization. The model with 
= 0.6, f^A = 0.4 is the output of the minimization code when made to fit the 
Q = 0.4, f^A = model. Temperature and polarization data were used for this fit. 
the two models differ in by 2. 

line (dashed) but not around Qm = 0.4. The squares and triangles correspond to the 
minima found by the minimization routine when constrained to move in subspaces of 
constant Q and agree with the ellipsoids from the Fisher matrix approach. 

Figure ^-9| shows the temperature and polarization spectra for the basis model and 
one found by the minimization routine with = (1495, .87, 0.6, 0.033, 0.051, 1.0, 0.39) 
where the last number is now Q^- This model differs from the basis model by a 

= 2 and so is practically indistinguishable from it. Only on large angular scales do 
the two models differ somewhat, but cosmic variance prevents an accurate separation 
between the two. In this case, polarization does not help to break the degeneracy. The 
agreement on the large angular scales is better for polarization than for temperature 
because the former does not have a contribution from the ISW effect, which is the 
only effect that can break this degeneracy. When both Q and ^a are included in the 
analysis the la error bars for both MAP and Planck increase. The greatest change 
is for the error bars on the curvature that now becomes AQ = 0.1 for both MAP 
and Planck. Note that improving the angular resolution does not help to break the 
degeneracy, which is why MAP and Planck results are similar. If one is willing to 
allow for both cosmological constant and curvature then there is a genuine degeneracy 
present in the microwave data and constraints from other cosmological probes will be 
needed to break this degeneracy. 
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5.4 Shape of the hkehhood function, priors and 
gravitational lensing 

As mentioned in Section 2, the Fisher information matrix approach used so far as- 
sumes that the hkehhood function is Gaussian around the maximum. In previous work 
0], this assumption was tested by calculating the likelihood along several directions 
in parameter space. This approach could miss potential problems in other directions, 
particularly when there are degeneracies between parameters. We will further test the 
Gaussian assumption by investigating the shape of the likelihood function varying one 
parameter at a time but marginalizing over the others. We fix the relevant parameter 
and let the minimization routine vary all the others in its search for the smallest x^- 
We then repeat the procedure for a different value of this parameter, mapping the 
shape of the likelihood function around the minimum. The minimization routine is 
exploring parameter space in all but one direction. These results may be compared 
with the prediction of the Fisher matrix which follow a parabola in the parameter 
versus log-likelihood plot. This comparison tests the gaussianity of the hkehhood 
function in one direction of parameter space at a time. 



The panels in Figure |5-10| show two examples of the results of this procedure. 
In most cases, the agreement between the Fisher matrix results and those of the 
minimization code is very good, especially very near the minimum (i.e., Ax^^ 2). 
As illustrated in the Qb panel, there are cases when increased more rapidly than 
predicted by the Fisher matrix. This is caused by the requirements that Uri, ^, 
and Qf, are all positive, which can be enforced easily in the minimization code. Of 
course, such priors are most relevant if the underlying model is very close to the 
boundary enforced by the prior and are only important on one side of the parameter 
space. The importance of this effect therefore depends on the underlying model. 
If the amount of information in the CMB data on a given parameter is sufficiently 
high, then the prior will have only a small effect near the maximum of the likelihood 
function. 

We also investigated the effect of gravitational lensing on the parameter recon- 
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Figure 5-10: Comparison between the Fisher matrix expansion of the hkehhood 
around the minimum (sohd hues) and direct minimization for two different cosmo- 
logical parameters. In most cases the agreement near the minimum is good. In the 
upper panel full triangles (crosses) correspond to fits of sCDM within the six param- 
eter family described in the text, including (not including) the effects of gravitational 
lensing. The lower curve belongs to the T/S = 0.28 model in the eight parameter 
space. In the lower panel, the > prior is reached for sCDM when Qf, > 0.05, 
which is why the minimization results differ from the Fisher matrix results. 



struction. As shown in |2^, gravitational lensing smears somewhat the acoustic 
oscillations but leaves the overall shape of the power spectrum unchanged. The am- 
plitude of the effect depends on the power spectrum of density fluctuations. Because 
the CMBFAST output consists of both CMB and density power spectra one can use 
them as an input for the calculation of the weak lensing effect following the method 
in 



. The gravitational lensing effect is treated self-consistently by normalizing the 
power spectrum for each model to COBE. We find that the addition of gravitational 
lensing to the calculation does not appreciably change the expected sensitivity to 
different parameters that will be attained with the future CMB experiments. This 
conclusion again depends somewhat on the underlying model, but even for sCDM 
where COBE normalization predicts two times larger small scale normalization than 
required by the cluster abundance data, the lensing effect is barely noticeable in the 
error contours for various parameters. 



5.5 Conclusions 

In this chapter, we have analyzed how accurately cosmological parameters can be 
extracted from the CMB measurements by two future satellite missions. Our work 
differs from previous studies on this subject in that we use a more accurate com- 
putational code for calculating the theoretical spectra and we include the additional 
information that is present in the polarization of the microwave background. We also 
investigate how the results change if we vary the number of parameters to be modeled 
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or the underlying model around which the parameters are estimated. Both of these 
variations can have a large effect on the claimed accuracies of certain parameters, so 
the numbers presented here should not be used as firm numbers but more as typical 
values. Of course, once the underlying model is revealed to us by the observations 
then these estimates can be made more accurate. The issue of variation of the errors 
on the number of parameters however remains, and results will always depend to 
some extent on the prior belief. If, for example, one believes that gravity waves are 
not generated in inflationary models (e.g. ||2^) or that they are related to the scalar 
perturbations through a simple relation (e.g. p4|), then the MAP errors on most 
parameters shrink by a factor of 2. Similarly, one may decide that models with both 
curvature and cosmological constant are not likely, which removes the only inherent 
degeneracy present in the CMB data. 

Using temperature data alone, MAP should be able to make accurate determina- 
tions (better than 10%) of the scalar amplitude (as), the baryon/photon ratio (fife/i^), 
the matter content (fi^/i^), the power spectrum slope {ug) and the angular diameter 
distance to the surface of last scattering (a combination of Q and Qa). If we restrict 
ourselves to models with no gravity wave content, then MAP should also be able 
to make accurate determinations of ^a, the Hubble constant and the optical depth, 
Kri- However, in more general models that include the gravity wave amplitude and 
spectral slope as additional parameters, the degeneracies between these parameters 
are large and they cannot be accurately determined. 

Several other measurements of the CMB anisotropics from the ground and from 
balloons are now in progress and accurate results are likely to be available by the 
time MAP flies. This additional information will help constrain the models further, 
especially determinations of the power spectrum at the smaller angular scales. 

Astronomical data can significantly reduce these degeneracies. The two nearly 
degenerate models, sCDM and a tilted vacuum dominated model (^a = 0.6) shown 
in Figure |5-3| can be distinguished already by current determinations of as^^'^, or by 
measurements of the shape of the galaxy power spectrum, or by measurements of the 
distance-magnitude relationship with SNela's. 
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map's measurements of polarization will significantly enhance its scientific return. 
These measurements will accurately determine the optical depth between the present 
and the surface of last scatter. This will not only probe star formation during the 
"dark ages" (5 < 2; < 1300), but will also enable accurate determination of the Hubble 
constant and help place interesting constraints on Qjy in models with tensor slope and 
amplitude as free parameters. 

Planck's higher sensitivity and smaller angular resolution will enable further im- 
provements in the parameter determination. Particularly noteworthy is its ability to 
constrain to better than 5% and the Hubble constant to better than 1% even in the 
most general model considered here. The proposed addition of polarization-sensitive 
bolometer channels to Planck significantly enhances its science return. Planck should 
be able to measure the ionization history of the early universe, thus studying pri- 
mordial star formation. Sensitive polarization measurements should enable Planck to 
determine the amplitude and slope of the gravity wave spectrum. This is particularly 
exciting as it directly tests the predicted tensor /scalar relations in the inflationary the- 
ory and is a probe of Planck scale physics. The primordial gravity wave contribution 
can at present only be measured through the CMB observations. One may therefore 
ask how well Planck can determine T/S assuming that other cosmological parameters 
are perfectly known by combining CMB and other astronomical data. The answer 
sensitively depends on reionization optical depth. Without reionization, T/S ~ 0.1 
can be detected, while with Uri = 0.1 this number drops down to T/S ^ 0.02. The 
equivalent number without polarization information is 0.2, regardless of optical depth. 
Improvements in sensitivity will further improve these numbers, particularly in the 
B polarization channel which is not cosmic variance limited in the sense that tensors 
cannot be confused with scalars. A detection of a. B component would mean a model 
independent detection of a stochastic background of gravitational waves or vector 
modes 

The most exciting science return from polarization measurements come from mea- 
surements at large angular scales (see Figures |5-6| and |5-7D . These measurements can 
only be made from satellites as systematic effects will swamp balloons and ground 
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based experiments on these scales. The low I measurements enable determinations 
of the optical depth and the ionization history of the universe and may lead to the 
detection of gravity waves from the early universe. Both foregrounds and systematic 
effects may swamp the weak polarization signal, even in space missions, thus it is 
important that the satellite experiment teams adopt scan strategies and frequency 
coverages that can minimize systematics and foregrounds at large angles. 

We explored the question of how priors such as positivity of certain parameters 
or constraints from other cosmological probes help reduce the uncertainties from the 
CMB data alone. For this purpose, we compared the predictions from the Fisher 
information matrix with those of the brute-force minimization which allows the easy 
incorporation of inequality priors. As expected, we find that positivity changes the er- 
ror estimates only on the parameters that are not well constrained by the CMB data. 
On the other hand, using some additional constraints such as the limits on the Hubble 
constant, age of the universe, dark matter power spectrum or qo measurements from 
type la supernovae can significantly reduce the error estimates because the degenera- 
cies present in these cosmological tests are typically different from those present in the 
CMB data. The minimization approach also allows testing the assumption that the 
log-likelihood is well described by a quadratic around the minimum, which is implicit 
in the Fisher matrix approach. We find that this is a good approximation close to the 
minimum, with no nearby secondary minima that could be confused with the global 
one. Finally, we also tested the effect of gravitational lensing on the reconstruction 
of parameters and found that its effect on the shape of the hkelihood function can be 
neglected. 

In summary, future CMB data will provide us with an unprecedented amount of 
information in the form of temperature and polarization power spectra. Provided 
that the true cosmological model belongs to the class of models studied here these 
data will enable us to constrain several combinations of cosmological parameters with 
an exquisite accuracy. While some degeneracies between the cosmological parameters 
do exist, and in principle do not allow some of them to be accurately determined 
individually, these can be removed by including other cosmological constraints. Some 
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of these degeneracies belong to contrived cosmological models, which may not survive 
when other considerations are included. The microwave background is at present our 
best hope for an accurate determination of classical cosmological parameters. 
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Table 5.1: Fisher matrix one-sigma error bars for different cosmological parameters 
when only temperature is included. Table 5.3 gives the cosmological parameters for 
each of the models. Columns with + correspond to MAP and those with x to Planck. 



Param. 


sCDM+ 


sCDM+ 


^ = 0.28+ 


^ =0.1+ 


^ =o.i>' 


Qa = 0.7+ 


Opcn+ 


Opcn*^ 




2.110-1 


4.210-1 


4.810-1 


4.710-1 


7.410-2 


4.110-1 


1.210-1 


4.710-2 


Ah 


1.710"^ 


9.210-^ 


1.110-1 


1.010-1 


5.1 10--' 


4.110-2 


2.010-2 


1.1 lo--' 


AHa 


9.810"^ 


5.310-1 


6.110-1 


5.810-1 


2.910-^ 


5.010-2 






AClih''' 


3.0 IQ-* 


1.010-^ 


9.8 10-* 


1.210-^ 


1.2 10-* 


9.7 10-* 


1.110-^ 


1.3 10-* 




1.210-1 


1.310-1 


1.410-1 


1.110-1 


8.210-^ 


1.910-1 


7.210-^ 


3.310-2 


Aus 


9.8 lO--' 


5.910-^ 


6.710-^ 


6.410-^ 


5.910-'' 


2.910-2 


2.410-^ 


3.310-'' 






3.910-1 


6.810-1 


5.310-1 


2.510-1 


3.210-1 






AriT 






3.910-1 


9.110-1 


9.410-1 


9.910-1 






An 














6.6 IQ-'^ 


5.2 10-* 







Table 5.2: Fisher matrix one-sigma error bars for different cosmological parameters 
when both temperature and polarization is included. Table 5.3 gives the cosmological 
parameters for each of the models. Columns with + correspond to MAP and those 
with X to Planck. 



Param. 


sCDM+ 


sCDM+ 


^ = 0.28+ 


^ =0.1+ 


^ =o.i>' 


Qa = 0.7+ 


Opcii+ 


Opcn^ 




1.8 10-2 


2.110-1 


2.810-1 


2.110-1 


1.010-2 


8.310-2 


6.510-2 


1.210-2 


Ah 


I.GIO - 


5.110 - 


5.8 10 - 


5.010 2 


3.010 


3.8 10 2 


1.9 10 2 


1.0 10-'' 


AQa 


9.310-^ 


2.910-1 


3.310-1 


2.910-1 


1.710-^ 


4.610-^ 






AClbh''' 


2.8 10-* 


6.1 10-* 


7.1 10-* 


6.2 10-* 


5.710-" 


8.9 10-* 


9.5 10-* 


1.1 10-* 




2.1 10-2 


2.110-2 


2.010-2 


2.010-2 


5.510-'' 


2.010-2 


3.210-2 


3.5 10-'' 


Aris 


4.810-^ 


3.110-2 


3.510-1 


3.010-2 


3.0 10-'' 


2.610-2 


1.710-2 


2.610-'' 






2.210-1 


4.310-1 


3.010-1 


4.510-2 


2.110-1 






Any 






3.910-1 


8.110-1 


1.710-1 


7.810-1 






AXe 










1.410-1 








An 














6.110-'^ 


4.1 10-* 







Table 5.3: Cosmological parameters for the models we studied. All models were 
normalized to COBE. 



Param. 


sCDM 


^ = 0.28 


W = 0-1 


Ha = 0.7 


Open 


h 


0.5 


0.5 


0.5 


.65 


.65 


Qa 


0.0 


0.0 


0.0 


0.7 


0.0 


n^ 


0.05 


0.05 


0.05 


0.06 


0.06 




0.05 


0.05 


0.1 


0.1 


0.05 


ris 


1.0 


0.96 


0.99 


1.0 


1.0 


T 


0.0 


0.28 


0.1 


0.0 


0.0 


riT 




0.04 


0.01 






n 


1.0 


1.0 


1.0 


1.0 


0.4 
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Appendix A 
Spin-weighted functions 



In this Appendix we review the theory of spin-weighted functions and their expansion 
in spin-s spherical harmonics. This was used in the main text to make an all-sky 
expansion of the Q and U Stokes parameters (Chapter ^). The main application of 
these functions in the past was in the theory of gravitational wave radiation (see e.g. 

Our discussion follows closely that of Goldberg et al. 0, which is based on the 
work by Newman and Penrose 0. We refer to these references for a more detailed 
discussion. 

For any given direction on the sphere specified by the angles [9, (p), one can define 
three orthogonal vectors, one radial and two tangential to the sphere. Let us denote 
the radial direction vector with n and the tangential with ci, 62. The latter two are 
only defined up to a rotation around n. 

A function sf{9,4') defined on the sphere is said to have spin-s if under a right- 
handed rotation of (61,62) by an angle ip it transforms as sf'{9,(f)) = e~*'^^ sf{9,(j)). 
For example, given an arbitrary vector a on the sphere the quantities a ■ 61 + ia ■ 62, 
n - a and a - ei — ia- 62 have spin 1,0 and —1 respectively. Note that we use a different 
convention for rotation than Goldberg et al. [0 to agree with the previous literature 
on CMB polarization. 

A scalar field on the sphere can be expanded in spherical harmonics, Yim{0,4'), 
which form a complete and orthonormal basis. These functions are not appropriate 
to expand spin weighted functions with s 7^ 0. There exist analogous sets of functions 
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that can be used to expand spin-s functions, the so called spin-s spherical harmonics 
sYim{0, (f)). These sets of functions (one set for each particular spin) satisfy the same 
completness and orthogonality relations, 



sY:^(e,<t>)sYim(e',<f>') = 5(<i>-<f>')5(cose-cose'). (a.i) 

Im 

An important property of spin-s functions is that there exists a spin raising (low- 
ering) operator d {d ) with the property of raising (lowering) the spin-weight of a 
function, {6 sf)' = e-'^'+^^'I'd sf, sf)' = e'^^"-^)'^^ J. Their explicit expression 
is given by 



^ sf{0,ct>) = -sin^(^) 
§ sfiOA) = -sin-^(e) 



(A.2) 



In this paper we are interested in polarization, which is a quantity of spin ±2. The 

d and d operators acting twice on a function ±2f{lJ^, <p) that satisfies 9^ sf = i'^ sf 
can be expressed as 



^"S/(/^,0) = (-9/.+ ^^) [(1 - //2) 2/(//» 



m 



:i-/i2) -2/(/^,0) 



(A.3) 



where ji — cos{9). With the aid of these operators one can express gYim in terms of 
the spin zero spherical harmonics YJ^, which are the usual spherical harmonics. 



Y' 



sYlm 




l^'Yim ,{0<S<1) 

-lyd-'Yi,^ ,{-i<s<o). 



(A.4) 
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The following properties of spin-weighted harmonics are also useful 



-y* 
Im 



s^lm 



[(/-S)(/ + S + 1)]^ 

-[{1 + S){l-S + 1)]'^ s-lYim 

-{I - S){1 + S + 1) sYlm 



(A.5) 



Finally, to construct a map of polarization one needs an explicit expression for the 
spin weighted functions, 



r 



sYim{n) = e 
'I - s 



im4> 



(/ + m)!(/-m)!2/ + lni/2 . 



r 



{l + s)\{l-s)\ Air 
' l + s -) 
r + s — ml 



sin'' (0/2) 



(A.6) 



For the special case |s| = 2 a more useful expression ( |2.13|) is provided in Chapter 
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